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ABSTRACT 


Measurement  and  analyses  of  the  thermal  response  of  five  identical  rocket 
nozzles  were  performed.  The  nozzles  were  fired  on  the  Philco-Ford  Aeronutronic 
Division,  TMSC-BM  simulator  between  April  15,  1971  and  May  4,  1971  under  Air 
Force  Contract  F04611-69-C-0039 .  The  simulated  solid  propellant,  designated 
AN3-3066 ,  contained  approximately  16  percent  aluminum.  The  five  nozzles  were 
exposed  to  various  multiple  pulse  duty  cycles  with  individual  pulse  times  rang¬ 
ing  from  10  to  20  seconds.  Nominal  chamber  pressure  for  each  pulse  was  750  psia 
Each  nozzle  was  instrumented  with  thermocoupled  plugs  at  two  locations  in  the 
exit  cone  and  with  two  spring  loaded  thermocouple  probes  at  the  backwall  of  the 
graphite  throat  package  parts. 

Post-test  measurements  of  surface  recession,  char  penetration,  density 
versus  depth,  and  in-depth  crystallization  were  made.  Predictions  of  the  ther¬ 
mal  and  ablative  responses  of  the  throat  package  and  the  e«it  cone  material  at 
an  area  ratio  of  approximately  4.0  were  performed  for  each  nozzle  and  pulse. 
Agreement  between  throat  package  backwall  temperature  response  data  and  the  pre¬ 
dicted  2-D  response  was  excellent.  One-dimensional  conduction  solutions  at  the 
throat  plane  were  less  accurate,  although  surface  temperature  and  recession 
predictions  agreed  favorably  with  the  2-D  conduction  results.  Recession  of  the 
PG  coated  throat  surface  was  underpredicted  using  the  kineticaily  controlled 
surface  response  model  developed  under  Air  Force  Contract  FCK 611-69-C-0031 .  The 
probable  reason  for  the  underpradiction  is  the  inability  to  accurately  evaluate 
the  severe  heating  at  the  boundary  layer  reattachment  point  or  the  throat  insert 

Temperature  response  predictions  for  the  FM  5055A  catbon  cloth  phenolic 
exit  cone  thermocoupled  plug  location  agreed  favorably  with  data,  although  con¬ 
duction  along  the  char  layer  during  nozzle  cooldowns  was  not  modeled.  Carbon 
phenolic  recession  comparisons  were  confounded  by  char  sv;ell  phenomenon.  Agree¬ 
ment  was  obtained  bv  utilizing  a  physically  reasonable  semi-empirical  char  swell 
correction  relatic  Density  profile  predictions  were  also  in  agreement  with 
X  ray  transmission  measurements  and  shave  and  weigh  measurements.  X-ray  dif¬ 
fraction  measurements  indicated  the  anticipated  increase  graphite  crystalliza¬ 
tion  with  charring,  independent  of  the  type  of  duty  cycle. 

The  general  conclusion  of  these  studies  is  that  current  thermal  analy¬ 
sis  techniques  are  suitable  for  the  evaluation  of  multiple  pulse  duty  cycle 
ablative  response  phenomena,  and  that  no  extraordinary  material  response  phe¬ 
nomena  occur  during  multiple  pulse  nozzle  exposures. 
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SECTION  1 


INTRODUCTION 

The  Solid  Rocket  Division  of  AFRPL  recently  conducted  an  in-house  program 
to  provide  design  information  for  multiple  pulse  duty  cycle  nozzles  in  solid 
rocket  motor  applications.  During  the  course  of  this  program,  five  instrumented 
nozzles  were  fabricated  and  tested.  The  material  components  for  these  nozzles 
included  a  MXC  313  carbon  phenolic  throat  entrance  section,  a  pyrolytic  graphite 
coated  nozzle  throat,  an  ATJ  graphite  throat  exit  section,  and  a  FM  5055A 
carbon  cloth  phenolic  exit  cone.  The  nozzles  were  tested  in  the  Aeronutronic 
Simulator  by  Philco-Ford  under  Contract  F04611-69-C-0039.  In  this  simulator, 
a  propellant  slurry  containing  approximately  16  percent  aluminum  was  used  to 
simulate  the  flame  temperature  and  combustion  products  of  ANB-3066  solid  propel¬ 
lant.  The  test  variables  included  in  the  multiple  pulse  duty  cycles  were  the 
number  of  pulses  which  varied  from  two  to  four,  the  length  of  each  pulse  which 
varied  from  10  to  20  seconds,  and  the  type  of  restart  which  included  both  hot 
and  cold.  The  nominal  chamber  pressure  for  each  pulse  was  approximately  750  psia. 

After  completing  the  test  firings,  the  nozzles  and  reduced  test  data 
were  delivered  to  Aerotherm  for  post-test  analysis  under  Contract  F04611-70-C-0019 . 
This  anlysis  included 

•  surface  recession  and  char  regression  measurements 

•  in-depth  density  profile  measurements  for  the  carbon  cloth  phenolic 
in  the  exit  cone  at  an  area  ratio  of  approximately  4.0 

•  in-depth  crystallographic  measurements  for  the  carbon  cloth  phenolic 
in  the  exit  cone  at  an  area  ratio  of  approximately  4.0 

•  prediction  of  transient  in-depth  temperatures  in  the  nozzle  throat 
insert,  a  comparison  of  predicted  and  measured  temperatures  at  the 
backside  of  the  AGSR  substrate  and  the  ATJ  exit  section  and  a  comparison 
of  predicted  and  measured  nozzle  throat  recessions 

•  comparison  of  predicted  and  measured  in-depth  transient  temperature 
distributions,  of  measured  and  predicted  char  regressions,  and  of 
measured  and  predicted  surface  recessions  for  the  carbon  cloth  phenolic 
located  in  the  exit  cone  at  an  area  ratio  of  approximately  4.0 

The  objectives  of  this  extensive  post-test  analysis  were  to  characterize 
the  degradation  of  the  charring  ablating  exit  cone  material  when  subjected  to 


multiple  pulse  duty  cycles,  to  obtain  material  thermal  performance  data  which 
can  be  effectively  used  by  the  nozzle  designer,  and  to  correlate  the  material 
thermal  performance  data  using  conventional  analytical  techniques.  The  analyti¬ 
cal  techniques  used  in  correlating  the  data  were  those  included  in  the  following 
computer  programs: 

•  Aerotherm  Chemical  Equilibrium  (ACE)  Computer  Program 

•  Aerotherm  Real  Gas  Energy  Integral  Boundary  (ARGEIBL)  Computer  Program 

•  Charring  Material  Ablation  (CMA)  Computer  Program 

•  Axisymmetric  Transient  Heating  and  Material  Ablation  (ASTHMA) 

Computer  Program 

These  analytical  techniques  were  chosen 

•  because  they  have  been  used  extensively  for  the  thermal  analysis  of 
rocket  nozzles  and  reentry  nose  cones 

•  because  they  have  been  validated  for  charring  ablating  and  carbon/ 
carbon  composite  materials  for  chamber  pressures  which  range  from 
200  psia  to  3000  psia  (References  1-5) 

•  because  they  are  currently  being  validated  for  pyrolytic  graphites 
under  Air  Force  Contract  F04611-69-C-0081 

•  because  they  are  or  shortly  will  become  available  to  the  rocket 
motor  community 

The  nozzle  configuration  used  in  the  multiple  pulse  duty  cycle  tests, 
the  test  facility  configuration,  and  a  summary  of  the  test  results  including 
chamber  pressure  and  temperature  histories  are  presented  in  Section  2.0.  A 
description  of  the  measured  material  performance  which  includes  the  surface 
recession,  char  regression,  in-depth  density  profiles,  and  x-ray  diffraction 
measurements  are  presented  in  Section  3.  A  description  of  the  analytical 
techniques  together  with  a  comparison  of  measured  and  predicted  material  thermal 
performance  is  presented  in  Section  4.  Finally,  conclusions  and  recommendations 
are  made  in  Section  5. 
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SECTION  2 


MULTIPLE  PULSE  DUTY  ROCKET  MOTOR  TESTS 

The  five  instrumented  nozzles  which  were  tested  in  the  Aeronutronic 
Solid  Rocket  Motor  Simulator  are  described  in  Section  2.1,  and  the  Simula'  r 
test  facility  is  described  in  Section  2.2.  The  propellant  slurry  used  to  simu¬ 
late  the  solid  rocket  motor  propellant  is  defined  in  Section  2.3.  The  duty 
cycles  including  the  chamber  pressure  histories  for  each  test  are  presented  in 
Section  2.4.  The  thermocouple  data  from  the  spring  loaded  thermocouples  behind 
the  nozzle  insert  and  from  the  plugs  in  the  nozzle  exit  cone  are  presented 
in  Section  2.5. 

2.1  NOZZLE 

The  five  nozzles  used  for  the  multiple  pulse  duty  cycle  and  the  geometry 
and  materials  of  each  nozzle  were  identical.  A  schematic  drawing  including 
the  average  nozzle  dimensions  is  shown  in  Figure  2-1,  and  the  detailed  pre-test 
quality  control  records  for  each  nozzle  are  summarized  in  Table  2-1.  As  shown 
in  Figure  2-1,  the  throat  entrance  consisted  of  MXC  313  carbon  phenolic  which 
overlaps  the  pyrolytic  graphite  (PG)  coating  at  the  entrance  to  the  nozzle  throat. 
The  PG  coating  which  has  a  nominal  thickness  of  0.049  mils  is  supported  by  an 
AGSR  graphite  substrate,  and  the  throat  insert,  is  insulated  from  the  steel 
shell  by  the  150  RPD  asbestos  phenolic.  The  throat  extensior  consists  of  ATJ 
graphite  backed  by  GA  carbon.  The  taper  on  the  backside  of  the  ATJ  graphite  to¬ 
gether  with  a  corresponding  taper  on  the  backside  of  the  GA  carbon  provided  the 
axial  structural  support  for  the  PG  coating  throat  insert.  The  exit  cone  is  com¬ 
posed  of  FM  5055A  carbon  cloth  phenolic  wrapped  parallel  to  the  nozzle  centerline. 

The  thermal  instrumentation  for  each  nozzle  consisted  of  two  spring  loaded 
Chromel/Alumel  thermocouples  and  two  thermocoupled  plugs.  The  pro-test  quality 
control  records  showed  that  all  of  the  thermocoupled  plugs  had  identical  dimensions 
and  that  the  locations  of  the  thermocouples  were  identical  for  all  of  the  plugs. 

The  dimensions  of  the  plug  and  the  location  of  the  three  in-depth  thermocouples 
are  shown  in  Figure  2-2.  As  shown  by  this  figure,  the  thermocouple  nearest  the 
surface  is  tungsten  5  percent  rhenium/tungsten  26  percent  rhenium  with  the  other 
two  being  Chromel/Alumel. 


TABLE  2-1 


PRE-TEST  NOZZLE  QUALITY  CONTROL  RECORD 


—  Dimensions  — 

(in) 

Ltr. 

Dim. 

Nozzle  Number* 

4 

5 

6 

7 

8 

Average 

A 

2.396 

2.402 

2.?02 

2 .400 

2.401 

2.400 

B 

3.S07 

3.509 

3.509 

3.504 

3.507 

3.507 

C 

4.754 

4.740 

4.755 

4.754 

4.755 

4.751 

D 

5.920 

5.926 

5.907 

5.910 

5.915 

5.915 

E 

7.336 

7.350 

7.336 

7.344 

7.338 

7.340 

F 

4.668 

4.662 

4.663 

4.663 

4.663 

4.663 

G 

2.294 

2.300 

2.304 

2.297 

2.296 

2.298 

H 

2.455 

2.457 

2.452 

2.455 

2.455 

kx: 

0.510 

0.520 

0.516 

0.521 

0.517 

mm 

0.050 

■ 

0.049 

0.051 

0.049 

U 

0.053 

B 

0.053 

0.055 

0.051 

The  nozzle  numbering  system  coincides  with  that  used  by  the 
fabricator  to  designate  the  nozzles  and  that  used  by  the  test 
facility  to  designate  the  tests. 


2 . 2  TEST  FACILITY 

The  tests  of  the  five  instrumented  nozzles  were  conducted  in  the  Aeronutronic 
EME  Solid  Propellant  Simulator  which  is  a  slurry  gaseous  rocket  engine  system. 

By  proper  selection  of  the  propellant  ingredients,  it  can  be  used  to  duplicate 
the  combustion  products  and  flame  temperature  of  aluminized  solid  rocket  propel¬ 
lants.  The  ingredients  in  the  simulation  propellant  for  the  multiple!  pulse  duty 
cycle  firings  consisted  of  oxygen,  hydrogen,  nitrogen,  and  a  slurry  mixture 
of  5  micron  particles  of  aluminum  and  aluminum  oxide  mixed  with  RP-1  fuel  and 
a  small  quantity  of  napalm.  These  propellant  ingredients  were  injected  into 
a  9.5  inch  diameter  combustion  chamber  which  was  lined  with  a  or.e-inch  thick 
hollow  cylinder  of  silica  phenolic.  Because  of  surface  recession  during  the  simu¬ 
lator  operation,  this  liner  was  replaced  after  conducting  tests  for  nozzles  4,  5, 
and  6.  The  surface  recession  on  the  new  liner  was  approximately  0.5  inches  at  the 
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conclusion  of  tests  for  nozzles  no.  7  anu  8.  Attached  to  the  simulator  combustion 
chamber  was  a  small  rocket  motor  which  was  used  to  assist  the  stcrtup  operation 
of  the  simulator.  The  start  sequence  was  as  follows: 

•  ignition  of  a  fuel  rich  mixture  of  oxygen  and  hydrogen  in  the  start 
motor  (The  combustion  products  were  injected  through  a  sonic  nozzle 
into  the  simulator  combustion  chamber.) 

•  after  approximately  0.75  seconds,  the  oxygen  and  hydcogen  propellant 
simulation  ingredients  were  injected  into  the  combu  .s:.ion  chamber 

•  after  approximately  0.25  seconds,  nitrogen  and  the  i. lurry  mixture 
were  also  injected  into  the  chamber 

•  after  approximately  1.0  second,  the  start  motor  was  shut  down 

During  each  test,  bcth  the  ballistic  and  thermcchemical  performance 
of  the  simulator  were  monitored.  The  ballistic  performance  was  determined  by 
pressure  measurements  made  in  the  simulator  combustion  chamber  and  in  the  start 
motor.  The  tnermochemical  performance  was  monitored  by  measuring  the  flowrates 
of  each  of  the  propellant  ingredients  which  were  then  compared  with  pre-selected 
values.  The  thermochemical  performance  of  the  simulator  in  terms  of  the  flame 
temperature  and  combustion  products  were  determined  from  theoretical  analysis 
using  the  pre-selected  flow  rates.  This  information  is  presented  in  the 
following  section. 

2.3  PROPELLANT 

As  mentioned  previously,  the  propellant  used  in  the  rocket  motor  simulator 
consisted  of  the  following  ingredients:  RP-1,  napalm,  five  micron  aluminum 
and  aluminum  oxide  particles,  oxygen,  hydrogen  and  nitrogen.  The  chemical  form¬ 
ula  and  nominal  relative  amounts  of  each  ingredient  are  shown  in  Table  2-2,  and 
the  resulting  chemical  elemental  concentrations  are  given  in  Table  2-3.  The 
flame  temperature  and  the  major  combustion  products  for  this  propellant  at  a 
chamber  pressure  of  700  psia  are  given  in  Table  2-4. 

The  propellant  formulation  given  in  Table  2-2  has  an  O/F  ratio  of  0.570 
where  the  O/F  ratio  is  defined  as 

o/F  =  - - -  (2-d 

^Ws lurry  +  Whydrogen) 

The  flowrates  and  O/F  ratio  for  each  test  are  defined  in  Table  2-5.  As  shown  by 
this  table,  the  O/F  ratio  varied  by  only  2  percent  (0.572  to  0.581). 
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TABLE  2-2 

SIMULATOR  PROPELLANT  INGREDIENTS 


Ingredient 

Chemical  Formula 

Relative 
Amount  by 
Mass 
Percent 

RP-l/Tri-Chlora- 

Ethylene 

C1.642H2.580°0.115N0.009C12.135 

30.4 

Napalm 

CH1. 9532 

5.13 

Aluminum 

A1 

11.31 

Aluminum  Oxide 

Al2°3 

8.24 

A11.0H35.P°5.0C13.0 

0.74 

Oxygen 

°2 

33.50 

Nitrogen 

N2 

8.70 

Hydrogen 

H2 

1.90 

TABLE  2-3 

SIMULATOR  PROPELLANT  ELEMENTAL  CONCENTRATIONS 


TABLE  2-4 


TABLE  2-5 


Another  parameter  often  used  to  characterize  a  solid  rocket  propellant 
is  the  oxidation  ratio  (XSO)  which  is  defined  by 


xs0  m  elemental  concentration  of  oxygen _ 

A  (elemental  concentration  of  oxygen  in  CO  +AL203) 

Applying  this  formula  to  the  simulator  propellant  gives  an  oxidation  ratio  of 
1.34  which  is  typical  of  solid  rocket  propellants  with  ar.  aluminum  loading  of  16 
percent. 

2.4  CHAMBER  PRESSURE  HISTORIES 

The  chamber  pressure  histories  for  the  five  multiple  pulse  duty  cycle 
rocket  nozzle  tests  are  described  in  Table  2-6  and  are  shown  graphically  in 
Figure  2-3.  More  exact  details  of  the  chamber  pressure  histories  are  given  in 
Appendix  A.  As  shown  by  this  figure,  the  number  of  pulses  per  test  varied  from 


TABLE  2-6 

DESCRIPTION  OP  MULTIPLE  PULSE  DUTY  CYCLES 


PRESSURE,  PSIA  PRESSURE,  PSIA  PRESSURE,  PSIA 


A*v»ne 


two  for  nozzle  no.  5  to  four  for  nozzle  no.  8,  the  pulse  duration  varied  from 
8  to  .10  seconds,  and  both  cold  and  hot  restarts  were  performed.  The  nominal 
chamber  pressure  for  all  the  pulses  remained  relatively  constant  as  it  varied 
only  from  680  to  750  psia.  In  addition,  the  total  test  time  was  approximately 
40  seconds  for  each  of  the  five  nozzles. 


2.5  MEASURED  IN-DEPTH  TEMPERATURES 

As  described  in  Section  2.1  and  shown  by  the  schematic  nozzle  drawing 
(refer  to  Figure  2-1) ,  the  nozzles  tested  during  this  program  had  1)  two 
spring  loaded  thermocouples  measuring  temperatures  at  the  backside  of  the  AGSR 
and  ATJ  graphites  in  the  nozzle  throat  assembly  and  2)  two  thermocoupled  plugs 
in  the  nozzle  exit  cone  with  three  thermocouples  in  each  plug.  The  maximum 
temperatures  measured  during  a  pulse  motor  burn  by  the  spring  loaded  thermocouples 
and  by  the  thermocouples  in  the  upstream  plug  are  summarized  in  Table  2-7.  De¬ 
tailed  transient  temperature  histories  for  all  of  the  thermocouples  are  presented 
in  Appendix  B. 


TABLE  2-7 

SUMMARY  OF  MEASURED  TEMPERATURES 


Measured  Temperatures  at  Motor  Burnout 

Max. 

Upstream  Thermocoupled  Plug  | 

Burn 

Backside 

Back side 

- ! 

Nozzle 

Pulse 

Tine 

of  AGSR 

of  ATJ 

TCI 

TC2 

TC3  I 

No. 

No. 

(sec) 

I'F) 

t’r) 

( *F) 

('F) 

CP)  , 

4 

1 

22 

900 

1500 

2200 

400 

150 

2 

18 

1000 

1330 

3000 

450 

100 

3 

9 

1880 

1840 

3260 

1200 

500 

5 

1 

22 

1070 

1300 

2480 

450 

150 

2 

22 

1950 

1600 

2400 

1250 

230 

6 

1 

17 

600 

1150 

2430 

320 

120 

2 

18 

1350 

ie2o 

3170 

1360 

520 

3 

11 

1830 

I860 

3240 

1670 

790 

7 

1 

17 

860 

1350 

1580 

250 

80 

2 

17 

1630 

1860 

2980 

860 

400 

3 

11 

2090 

2180 

3120 

1500 

700 

8 

1 

12 

630 

1230 

930 

150 

— 

2 

12 

1150 

1570 

2850 

1350 

3 

12 

1620 

2030 

2380 

650 

— - 

4 

13 

2090 

2300 

3070 

1860 

..  _ 

By  comparing  measured  temperatures  presented  in  Table  2-7  for  similar 
pulse  durations  and  for  identical  restart  conditions,  an  estimate  can  be  made  of 
the  accuracy  of  the  temperature  measurement.  Similar  environmental  conditions 
exist  at  the  surface  of  the  nozzle  for  pulse  no.  1  for  nozzles  no.  4  and  5  and 
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and  for  pulse  no.  1  for  nozzles  no.  6  and  7.  For  nozzles  no.  4  and  5,  pulse 
no.  1  has  a  duration  of  22  seconds  while  for  nozzles  no.  6  and  7,  pulse  no.  1 
has  a  duration  of  17  seconds.  For  the  first  two  nozzles ,  the  measured  tempera¬ 
tures  are  900  and  1070°F  at  the  backside  of  the  AGSR  and  1500°F  at  the  backside 
of  the  ATJ.  These  variations  in  the  measured  temperatures  give  an  accuracy  of 
approximately  ±  100°F  for  the  spring  loaded  C/A  thermocouples.  Making  a  similar 
comparison  for  the  upstream  thermocoupled  plug  shows  that  essentially  identical 
temperatures  are  measured  for  pulse  no.  1  for  nozzles  no.  4  and  5.  However,  a 
large  discrepency  exists  in  the  measured  temperatures  for  nozzles  no.  6  and  7. 
This  discrepancy  which  is  approximately  850°F  for  TC  no.  1  is  probably  a  result 
of  an  error  in  defining  the  thermocouple  depth.  Based  on  this  comparison, 
the  accuracy  of  the  temperature  measurement  is  approximately  +  400 °F  for  TC  no.  1 
Because  of  the  lower  temperature  gradients,  the  magnitude  of  the  temperature 
variation  would  be  less  at  the  other  two  thermocouple  locations  in  the  thermo¬ 
coupled  plug. 
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SECTION  3 


POST-TEST  MEASURED  MATERIAL  PERFORMANCE 

The  post-test  analyses  to  define  the  material  thermal  response  of  the 
multiple  pulse  duty  cycle  nozzles  included: 

•  measurement  of  surface  recession  and  char  regression  profiles 
(Section  3.1) 

•  measurement  of  in-depth  material  density  profiles  for  the  FM  5055A 
carbon  cloth  phenolic  at  a  supersonic  area  ratio  of  approximately  4.0 
(Section  3.2) 

•  measurements  of  in-depth  crystal  structure  profiles  for  the  FM  5055A 
carbon  cloth  phenolic  at  a  supersonic  area  ratio  of  approximately  4.0 
(Section  3.3) 

3.1  SURFACE  RECESSION  AND  CHAR  THICKNESS  MEASUREMENTS 

Surface  recession  and  char  thickness  data  were  obtained  at  several  loca¬ 
tions  along  each  of  eight  circumferential  planes  in  the  five  multiple  pulse 
nozzles.  The  measurement  planes  were  spaced  evenly  at  every  45°  of  rotation 
about  the  nozzle  centerline  with  the  0  degree  reference  plane  passing  through 
the  thermocoupled  plugs.  Diameter  measurements  were  taken  by  AFRPL  personnel  at 
one  circumferential  location  but  at  several  nozzle  axial  locations  to  confirm 
the  pre-test  surface  geometry.  In  defining  the  post  test  surface  recession, 
the  assumption  was  made  that  the  initial  internal  nozzle  shape  was  axisymmetric 
and  that  the  initial  radius  was  half  the  measured  diameter  in  all  circumferential 
planes.  The  surface  recessions  were  determined  by  subtracting  the  post-test 
measurements  of  surface  radius  from  the  corresponding  pre-test  dimension. 

The  post-test  radii  were  defined  using  the  technique  described  in 
Reference  6.  The  technique  consisted  of  chucking  the  entire  nozzle  assembly 
in  a  machine  shop  lathe  so  that  the  lathe  tool  position  indicator  could  be  used 
to  define  the  distance  from  the  lathe  (and  nozzle)  centerline  to  the  surface 
point  of  interest  (i.e.  the  local  radius).  Measurements  of  nozzle  radius  at 
several  discrete  surface  points  in  one  plane  defined  the  post-test  surface 
contours  and  enabled  the  evaluation  of  the  local  surface  recession.  The  uncer¬ 
tainty  in  radius  measurements  taken  in  this  manner  was  estimated  to  be  +  0.002 
inches.  The  post-test  nozzle  surface  contours  determined  using  this  measurement 
technique  are  compared  to  pre-test  contours  in  Figures  3-1  through  3-5.  Photo¬ 
graphs  showing  typical  surface  characteristics  are  shown  ir  Figure  3-6. 


LINE  DESCRIPTIONS 
ON  FIGURE  3-la. 


f)  225°  PLANE 
FIGURE  3-1  (CONTINUED) 


LINE  DESCRIPTIONS 
ON  FIGURE  3-2a . 


CHAR  LINE 


a)  0°  PLANE  (THERM0C0UPLED  PLUG  PLANE) 


FIGURE  3-3  PRE,  POST,  AND  CHAR  LINE  CONTOURS  FOR  NOZZLE  NO.  6 


b)  45°  PLANE 

FIGURE  3-4  PRE,  POST,  AND  CHAR  LINE  CONTOURS  FOR  NOZZLE  NO.  7 


PRE  TEST  CONTOURS 


/ 

/ 


POST 


TEST  CONTOURS 


a)  0°  PLANE  (THERMOCOUPLED  PLUG  PLANE) 


b)  45°  PLANE 


FIGURE  3-5  PRE,  POST,  AND  CHAR  LINE  CONTOURS  FOR  NOZZLE  NO.  8 


d)  135°  PLANE 


FIGURE  3-5  (CONTINUED) 
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Figure  3- 


j)  i;xit  Cone 


Section 


Typical  Multiple  Pulse  Nozzle  Post  Test  Surface 
Characteristics,  Nozzle  No.  8. 


SIPS 


Pre-  and  post-test  radii  and  net  surface  recession  at  the  two  axial  loca¬ 
tions  (nozzle  throat  and  upstream  thermocouple  plug)  for  which  analysis  efforts 
were  performed  are  given  in  Tables  3-1  and  3-2.  The  throat  recession  data  pre¬ 
sented  in  Table  3-1  are  plotted  as  a  function  of  circumferential  location  in 
Figure  3-7.  Also  included  in  Table  3-1  are  average  throat  recession  rate  values 
per  pulse  which  have  been  inferred  from  chamber  pressure  versus  time  data. 

Both  the  surface  contours  and  the  nozzle  radius  tabulations  show  that 
for  all  nozzles,  except  nozzle  number  5,  the  carbon  cloth  phenolic  in  the  exit 
cone  exhibited  some  amount  of  surface  growth  (negative  recession) .  To  verify 
that  these  data  were  not  grossly  in  error  due  to  inaccurate  post- test  radius 
measurements,  the  recessions  were  evaluated  by  an  alternate  procedure.  After 
the  nozzle  radius  measurements  were  completed,  the  ablative  components  were 
separated  from  the  steel  shells  and  sectioned  along  the  eight  circumferential 
planes.  Measurements  of  the  final  exit  cone  material  thickness  at  several  loca¬ 
tions  in  nozzles  no.  5  and  8  were  made,  and  net  surface  movements  were  computed. 
The  results  obtained  using  these  measurements  were  within  0.010  inches  of  those 
obtained  from  the  lathe  data.  Possible  explanations. for  the  surface  growth 
phenomena  are  presented  in  Section  4. 


Air.u  nx-  rutfixocourito  nur.  plane,  cegp.ces 

FIGURE  2-7  CORRELATION  OF  LOCAL  NOZZLE  THROAT  SURFACE  RECESSION 
WITH  CIRCUMFERENTIAL  LOCATION 
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TAILS  3-1 


SUMMARY  Of  SOULS  THROAT  RECESSION  HMSURSNMTS 


Resale 

Rusher 

_ 

hhp 

mm 

!HH1 

Aveng. 

R*o.(«lcoc 

Firing 

M— 

K«c«asio 

Inftrrtd 

Fx«i*ur« 

(Bili/f 

0** 

45* 

90* 

135* 

315* 

m 

Rost -Teat  Radius,  laches 

l.sio 

1.255 

1.245 

B9 

1.145 

1.140 

1.245 

KH 

i 

2.0 

rrs-Tsst  Radius,  inches 

1.147 

1.147 

1.147 

1.147 

1.147 

1.147 

Intll 

2 

2.1 

Surtees  Rsoassion 

(#343) 

(.104) 

(.0*1) 

(.0*3) 

am 

-.007 

(.0*8) 

(.151) 

.037  ft  .050 

3 

3.7 

-InH 

1.170 

1.170 

1.170 

1.170 

Iwrf 

1.142 

1.170 

1.144 

i 

l.» 

■i 

1.150 

1.150 

1.150 

1.150 

jgBl 

1.150 

1.150 

1.150 

2 

0.4 

;  j|i|l§ 

.020 

.020 

.020 

.020 

.005 

.012 

MpTfll 

.012 

jj$9|§ 

1.315 

1.395 

1.140 

1.155 

1.150 

IWTfH 

1.2(5 

i 

1.4 

1.152 

1.152 

1.152 

1.152 

1.152 

1.152 

2 

2.C 

1 

(.233) 

(.233) 

.020 

-.002 

ieui 

(.143) 

(.113) 

.035  t  .030 

3 

0.9 

1.445 

1.190 

1.200 

IW- 1 

1.170 

1.1(5 

1.345 

1 

2.3 

ifSiMn 

1.143 

1.143 

1.143 

l»«l 

Hi 

1.143 

1.143 

1.143 

2 

2.3 

(.303) 

(.057) 

.017 

.027 

.042 
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Char  thickness  measurements  were  taken  from  the  sectioned  pieces  of  the 
ablative  exit  cone  with  the  char  thickness  defined  as  the  distance  from  the  in¬ 
ternal  surface  to  the  midpoint  of  the  char/virgin  material  interface.  Uncertainty 
in  this  measurement  is  due  primarily  to  the  difficulty  in  visually  picking  out 
the  interface  since  both  the  virgin  material  and  the  char  are  black.  The  uncer¬ 
tainty  in  these  data  points  is  estimated  to  be  +  0.020  inches.  Sufficient  char 
thickness  measurements  were  acquired  to  define  a  char  line  contour  for  each  cir¬ 
cumferential  plane  of  the  five  nozzles.  These  contours  are  shown  along  with  pre- 
and  post-test  surface  contours  in  Figures  3-1  through  3-5?  in  addition,  the  thick¬ 
ness  measurements  at  the  A/A*  >4.2  prediction  location  are  given  in  Table  3-1. 

3.2  IN-DEPTH  DENSITY  PROFILE  MEASUREMENTS 

The  in-depth  material  density  profiles  were  measured  for  the  FM  5055a 
carbon  cloth  phenolic  at  a  supersonic  area  ratio  of  4.2.  This  area  ratio  cor¬ 
responds  to  the  location  of  the  upstream  thermocoupled  plug.  The  two  techniques 
used  for  measuring  the  density  profiles  were  the  shave  and  weigh  technique  and 
the  x-ray  transmission  technique.  For  the  shave  and  weigh  determination  of 
density  profiles,  two  samples  were  machined  from  the  exit  cone  material  of  each 
nozzle.  Samples  were  taken  from  the  0  and  180  degree  circumferential  locations 
at  the  area  ratio  of  4.2.  Each  sample  was  0.5  inches  square  in  the  plane  cor¬ 
responding  to  the  exposed  surface  and  extended  approximately  0.75  inches  in  the 
direction  normal  to  the  surface  (i.e.  to  the  steel  shell).  Density  as  a  func¬ 
tion  of  depth  was  evaluated  by  a  differencing  technique.  Overall  material  den¬ 
sity  was  accurately  obtained  from  dimension  and  weight  measurements.  Then,  ap¬ 
proximately  25  mils  of  the  sample  was  removed  starting  from  the  fully  charred, 

0.5  inch  square  end?  and  the  sample  was  re-measured.  The  change  in  weight  and 
volume  defined  the  density  of  the  removed  material.  This  process  was  repeated 
until  the  sample  was  consumed.  The  density  profiles  resulting  from  these 
measurements  are  represented  by  the  symbols  in  Figure  3-8.  The  curve  shown  on 
each  of  these  figures  is  the  analytical  prediction,  and  it,  as  well  as  the  com- 

Par**on  between  the  measured  and  predicted  profiles,  is  discussed  in  Section 
4.3.3. 

The  second  method  used  to  define  the  density  profile  was  the  x-ray  trans¬ 
mission  technique.  The  experimental  portion  of  this  technique  is  performed  by 
directing  a  narrow  x-ray  beam  of  known  intensity  normal  to  one  surface  of  the 
material  sample  and  detecting  the  intensity  of  the  beam  which  is  transmitted 
through  the  material  and  out  the  opposite  surface.  The  configuration  of  this 
apparatus  is  shown  schematically  in  Figure  3-9.  The  x-ray  source  generates  a 
1  mm  wide  beam  of  ka  radiation  which  is  passed  through  a  kB  nickel  filter  and 
then  through  a  collimator  which  directs  it  at  the  sample  being  tested.  The 
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portion  of  the  beam  which  passes  through  the  sample  is  measured  by  the  detector 
and  recorded  on  a  strip  chart  recorder.  The  sample  was  positioned  in  the  x-ray 
beam  at  discrete  1  mm  steps.  The  procedure  for  making  the  measurement  was  to 
1)  locate  the  sample,  2)  turn  on  the  strip  chart  recorder  which  had  a  speed 
of  1/2  inch/minute,  3)  turn  on  the  x-ray  source,  4)  record  data  for  2  minutes. 
The  x-ray  source  and  recorder  were  then  stopped  and  the  sample  was  repositioned. 
Measurements  were  taken  at  1  mm  intervals  which  also  corresponded  to  the  diameter 
of  the  x-ray  source. 

The  material  sample  was  taken  from  nozzle  number  8  at  the  location  shown 
in  Figure  3-10.  The  nominal  dimensions  of  the  sample  are  also  shown. 


FIGURE  3-10  LOCATION  FOR  X-RAY  TRANSMISSION  SAMPLE  (NOZZLE  NO.  8) 


The  relative  transmission  measurements  of  the  x-ray  beam  through  this  sample 
are  shown  in  Figure  3-11  for  the  various  discrete  locations.  The  abscissa  on  this 
figure  is  distance  from  the  surface  which  corresponds  to  the  center  of  the  x-ray 
source;  thus,  at  a  distance  corresponding  to  zero,  half  of  the  beam  is  passing 
through  air  and  half  through  the  sample.  The  data  shown  in  Figure  3-11  were 
reduced  to  material  density  using  the  following  equation: 


j  =  exp 
o 

where 

I  -  the  transmitted  x-ray  intensity 
IQ  -  the  x-ray  incident  beam  intensity 


(3-1) 


^  -  the  mass  absorption  coefficient 

p  -  local  material  density 

x  -  the  thickness  of  the  absorbing  material 
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DISTANCE  FROM  SURFACE,  MM 
FIGURE  3-11  MEASURED  X-RAY  TRANSMISSION 
FOR  DENSITY  PROFILE  SAMPLE 


A  mass  absorption  coefficient  of  1.15  cm2/gr  was  evaluated  based  on  the  transmis¬ 
sion  intensity  through  the  virgin  material.  The  virgin  material  density  used 
for  this  evaluation  was  1.37  gr/cm3.  This  value  was  obtained  from  a  small  sample 
of  virgin  material  which  was  taken  from  tne  nozzle  in  the  vicinity  of  the  x-ray 
sample.  The  mass  absorption  coefficient  for  the  virgin  material  carbon  cloth 
phenolic  compares  favorably  with  the  value  (1.18  cm2/gr)  obtained  from  a  sample 
of  pure  graphite  for  the  x-ray  transmission  results  presented  in  Reference  1. 

The  densities  calculated  using  the  above  equation  are  shown  in  Figure  3-12 
The  amount  of  scatter  in  the  data  is  a  result  of  the  porous  nature  of  the  sample 
which  is  shown  by  the  photograph  presented  in  Figure  3-13.  In  addition,  the 
porosity  in  the  char  region  probably  explains  the  low  value  of  char  density  deter¬ 
mined  from  the  x-ray  measurements  in  the  char  zone. 

3.3  X-RAY  DIFFRACTION  MEASUREMENTS 

X-ray  diffraction  measurements  were  made  for  three  samples  of  carbon  cloth 
phenolic  in  each  of  the  five  nozzles.  As  shown  in  Figure  14  th  •  three  samples 
in  each  of  the  nozzles  were  one-eighth  inch  thick  and  one-half  inch  on  a  side 
and  were  taken  at  an  area  ratio  of  4.2  and  at  the  180  degree  circumferential 
location.  The  three  samples  were  machined  from  the  exit  cone  at  depths  corres¬ 
ponding  to  the  fully  charred  zone,  the  partially  charred  zone,  and  the  virgin 
material  zone. 
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The  reason  for  making  the  x-ray  diffraction  measurements  was  to  determine 
whether  any  significant  changes  occurred  in  the  crystal  structure  of  the  material 
as  a  result  of  subjecting  the  material  to  multiple  pulse  duty  cycles.  The  type 
of  material  changes  which  would  be  expected  is  an  increase  in  the  graphitic  carbon 
content  in  the  char.  This  increase  in  graphitic  carbon  can  be  measured  by  x-ray 
diffraction  oecause  graphitic  carbon  {or  graphite)  has  a  uniformity  of  structure 
which  can  be  represented  by  the  graphite  crystal  lattice  shown  in  Figure  3-15. 


As  this  type  of  crystal  structure  becomes  more  prevalent,  the  x-ray  diffraction 
measurement  will  begin  to  show  a  sharp  peak  at  an  incidence  angle  defined  by 
Bragg's  Law 


sin  O  =  ^ 


(3-2) 
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where 


0  -  angle  of  incidence 

X  -  wavelength  of  x-ray  source 

d  -  spacing  between  adjacent  crystal  planes 

For  the  x-ray  diffraction  measurements  made  under  this  program,  a  copper  ka 

0 

x-ray  source  was  used  which  has  a  wavelength  of  1.54  A.  When  this  value  for 

O 

lambda  is  used  in  the  above  equation  together  with  a  value  of  3.355  A  for  the 
spacing  between  adjacent  planes  of  carbon  atoms  in  the  graphite  crystal  lattice, 
the  angle  of  x-ray  incidence  which  gives  the  maximum  strength  of  the  diffracted 
x-ray  beam  is  13.3  degrees. 

The  results  from  the  x-ray  diffraction  measurements  for  the  nozzle  material 
samples  are  shown  in  Figure  3-16.  The  ordinate  for  each  of  these  figures  repre¬ 
sents  on  a  relative  basis  the  strength  of  the  diffracted  x-ray  beam.  The  abscissa 
represents  the  angle  20  (two  times  the  incidence  angle  defined  in  Equation  (3-2)). 
The  samples  for  nozzle  no.  8  (Figure  3-16e)  were  tested  first.  The  incident 
angle  of  the  x-ray  beam  was  varied  from  45  (or  26  =  90°)  to  5  degrees  to  determine 
the  range  of  interest  for  this  angle.  As  shown  by  this  figure,  the  x-ray 
diffraction  measurements  were  of  primary  interest  for  2Q  angles  between  10  and 
40°,  and,  therefore,  sample  sets  for  the  other  nozzles  were  only  tested  over  this 
range.  For  all  of  the  sample  sets,  the  x-ray  diffraction  measurements  for  the 
virgin  material  sample  gave  a  low  magnitude  increase  over  a  broad  angle  spectrum 
from  between  10  and  30  degrees.  For  the  partially  charred  .samples,  very  little 
change,  except  for  nozzle  no.  5,  was  noted  between  the  virgin  and  partially 
charred  x-ray  diffraction  measurements.  However,  for  the  fully  charred  sample, 
the  intensity  of  the  diffracted  beam  increased  for  20  angles  between  22  and  28 
degrees  with  the  peak  occurring  at  a  20  incident  angle  of  approximately  26  degrees 
These  changes  in  the  x-ray  diffraction  measurements  for  the  fully  charred  samples 
indicate  that  crystallization  and  ordering  of  the  (00.2)  layer  plane  of  the 
hexagonal  graphite  crystal  has  taken  place.  The  most  pronounced  peaks  occurred 
for  samples  from  nozzles  No.  4  and  6  (Figures  3-16a  and  3-16c) . 

Since  the  amount  of  graphitic  carbon  content  is  a  function  of  the  maximum 
temperature  and  of  the  time  at  this  temperature,  the  trends  shown  by  the  x-ray 
diffraction  measurements  for  nozzles  no.  4  and  6  were  compared  with  the  predicted 
surface  temperatures  and  with  the  temperatures  measured  by  TC  no.  1.  This  compari 
son  showed  that 
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e)  NOZZLE  NO.  8 
FIGURE  3-16  (Concluded) 

•  the  maximum  temperatures  recorded  by  TC  no.  1  were  for  nozzles 

no.  4  and  6  (Table  2-7)  which  is  consistant  with  the  x-ray  diffraction 
measurements 

•  the  maximum  predicted  surface  temperatures  and  the  duration  at 
this  temperature  are  not  significantly  different  for  the  five  duty 
cycles  (Figures  4-25  through  4-29)  which  is  inconsistant  with  the  x-ray 
diffraction  measurements 

Because  the  magnitude  of  the  measured  temperatures  in  the  thermocoupled  plug  is 
a  strong  function  of  the  thermocouple  depth  and  because  small  variations  in  this 
depth  would  result  in  significant  changes  in  the  measured  temperature,  the 
apparent  correlation  between  the  x-ray  diffraction  and  the  thermocouple  measure¬ 
ments  should  be  viewed  with  caution.  This  is  particularly  true  because  of  the 
discrepancies  in  the  measured  temperatures  noted  in  Section  2.5  and  because 
the  predicted  temperatures  do  not  substantiate  this  same  correlation. 
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EVALUATION  OF  MATERIAL  PERFORMANCE  DATA 

The  simulator  firings  of  the  five  nozzles  provide  information  for  the 
design  of  multiple  pulse  duty  cycle  rocket  motors.  In  this  section,  the  multiple 
pulse  firing  data  are  used  to  verify  the  adequacy  of  current  analysis  procedures. 
Once  substantiated,  these  analysis  tools  should  facilitate  a  more  efficient  design 
design  process  by  minimizing  the  amount  of  experimentation  required.  Also  in 
this  section  the  material  response  predictions  are  used  to  identify  the  relative 
significance  of  observed  trends  in  the  firing  data.  Any  material  response 
characteristics  peculiar  to  multiple  pulse  duty  cycles  would  be  exposed  by  the 
correlation  effort.  Results  of  response  analyses  are  compared  to  material  per¬ 
formance  data  for  each  of  the  five  multiple  pulse  nozzles. 

The  methods  utilized  to  perform  these  analyses  are  summarized  in  Section 
4.1.  The  important  boundary  layer,  surface  thermochemistry,  and  material  prop¬ 
erty  information  which  are  needed  for  the  material  response  calculations  are 
described  in  Section  4.2.  Results  of  the  response  analyses  are  compared  to  the 
test  data  in  Section  4.3. 

4.1  SUMMARY  OF  ANALYSIS  METHODS 

The  analysis  of  material  response  to  rocket  nozzle  environment  encompasses 
both  the  determination  of  surface  recession  due  to  chemical  corrosion  and  mechan¬ 
ical  erosion  and  the  evaluation  of  in-depth  heat  conduction  and  material  decom¬ 
position.  The  events  occurring  at  and  below  the  surface  of  a  charring  ablator 
which  is  being  exposed  to  a  solid  propellant  rocket  environment  are  characterized 
in  Figure  4-1.  Except  for  the  material  decomposition  and  pyrolysis  gas  flow,  the 
figure  applies  equally  to  noncharring  materials  such  as  pyrolytic  graphite.  A 
set  of  comprehensive  computerized  analysis  tools  have  been  developed  to  define 
the  thermochemical  performance  of  these  types  of  materials.  These  computer  codes 
are  used  to  evaluate  thermal  ablation  phenomena,  and  they  have  been  designed  to 
be  applicable  for  arbitrary  material  and  propellant  compositions  and  for  any  gen¬ 
eral  environmental  duty  cycle  specification. 
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The  ablation  analysis  computer  codes  utilized  for  the  multiple  pulse 
nozzle  calculations  are  as  follows: 

•  Aerotherm  Chemical  Equilibrium  (ACE)  Program* 

•  Aerotherm  Real  Gas  Energy  integral  Boundary  Layer  (ARGEIBL)  Program 

•  Charring  Material  Ablation  (CMA)  Program 

•  Axi-Symmetric  Transient  Heating  and  Material  Ablation  (ASTHMA) 
Program 

The  first  program  computes  the  thermochemical  state  of  any  set  of  chemical  ele¬ 
mental  quantities.  For  ablative  rocket  nozzle  considerations ,  the  code  assists 


The  ACE  program  combines  all  the  features  of  the  Equilibrium  Surface  Thermo¬ 
chemistry,  Version  3  (EST3)  code  (currently  available  to  the  rocket  technology 
community)  and  the  graphite  response  code  developed  under  Contract  F04611-69-C 
0081  to  account  for  surface  kinetics. 
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in  the  computations  of  nozzle  expansion  properties,  gas  phase  transport  and 
thermodynamic  state  properties,  pyrolysis  gas  enthalpy-temperature  relation¬ 
ships,  and  thermochemical  behavior  of  a  surface  material  in  contact  with  the 
hyperthermal  boundary  layer  flow.  The  surface  removal  mechanisms  considered 
by  the  code  are  chemical  corrosion,  surface  vaporization,  and  liquid  layer  re¬ 
moval. 

The  ARGEIBL  program  calculates  bulk  film  transfer  coefficients  for  non¬ 
ablating  laminar  and/or  turbulent  boundary  layers  of  axisymmetric  or  flat  plate 
geometries  using  an  energy  integral  method.  The  solution  procedure  allows  any 
real  gas  chemical  system  because  it  accepts  generalized  thermodynamic  property 
information  which  are  obtained  from  ACE  code  computations.  Variations  in  axial 
surface  temperature  and  free  stream  properties  are  accounted  for  in  the  integral 
procedure.  The  convective  heat  transfer  coefficient  distributions  computed  us¬ 
ing  the  ARGEIBL  technique  have  been  confirmed  by  solutions  generated  with  the 
more  sophisticated  Boundary  Layer  Integral  Matrix  Procedure  (BLIMP)  computer 
code  (Ref.  5) . 

Both  the  ACE  and  ARGEIBL  codes  are  utilized  to  provide  information  needed 
by  the  CMA  or  the  ASTHMA  material  response  codes.  The  CMA  code  computes  the 
transient  thermal  response  of  a  material  which  is  receding  due  to  chemical  or 
mechanical  surface  removal  and  which  is  decomposing  in  depth  via  specified  rate 
controlled  phase  change  reactions.  Heat  conduction  is  modeled  along  a  one¬ 
dimensional  nodal  column  of  arbitrarily  varying  cross-sectional  area  and  is  im¬ 
plicitly  tied  to  ablating  or  nonablating  surface  boundary  conditions.  The  ASTHMA 
code  is  used  to  evaluate  the  transient  response  of  nondecomposing  materials  such 
as  graphite  within  which  two-dimensional  heat  conduction  effects  are  important. 
The  same  generalized  surface  boundary  condition  options  used  by  CMA  are  avail¬ 
able  in  ASTHMA . 

The  major  features  of  the  four  computer  codes  are  described  further  in 
Table  4-1  and  the  procedure  for  using  these  codes  to  perform  an  ablation  and 
thermal  response  analysis  is  diagrammed  in  Figure  4-2.  The  sequence  of  events 
in  the  analysis  procedure  is  also  outlined  in  Table  4-2.  The  major  result  of 
each  analysis  step  in  the  procedure  is  defined  along  with  how  this  result  is 
used  to  obtain  the  desired  analysis  predictions.  The  analysis  starts  by  defin¬ 
ing  the  boundary  layer  edge  properties  from  the  flow  field  and  ACE  expansion 
analyses  and  concludes  with  the  definition  of  the  thermal  performance  parameters, 
namely,  char  depth,  surface  recession,  and  in-depth  temperature  distributions. 
Additional  details  of  the  computerized  analysis  techniques  utilized  in  the  pro¬ 
gram  are  included  in  Appendix  c. 
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TABLE  4-1 


SUMMARY  Or  MAJOR  FEATURES  OF  COMPUTER  PROGRAMS 


Computer 

Program 

Major  Features  j 

Typical  Use 

|  ACS  "”~ 

1. 

Calculation  of  thermodynamic  state  for  ! 
isentropic  or  isoenthalpic  process  for 
a  closed  system  including  condensed 
phases . 

Kozsle  equilibrium  isentropic 
expansion. 

1 

1 

i 

2. 

Calculation  of  transport  properties 
(viscosity,  thermal  conductivity,  and 
mass  diffusion  coefficients). 

Evaluation  of  transport  coeffi¬ 
cients  for  use  in  boundary 
layer  analysis. 

1 

i 

1 

i 

3. 

Calculation  of  surface  state  of  any 
material  exposed  to  an  environment  con¬ 
sidering  all  possible  chemical  reac¬ 
tions  between  the  material  and  the 
environment.  Surface  reactions  may  be 
kinetically  controlled.  Mechanical 
removal  of  surface  melt  layer  may  also 
be  considered. 

Generate  Mollier  charts  for 
any  gas  system. 

• 

4. 

Considers  unequal  diffusion  coeffi¬ 
cients. 

Generation  of  thermochemical 
data  for  use  with  CMA  and 

ASTHMA  computer  programs. 

;  CMA 

i 

i 

1. 

Calculation  of  transient  thermal  re¬ 
sponse  of  a  composite  slab  or  cylinder 
containing  up  to  five  in-depth  charring 
(or  noncharring)  materials. 

Thermal  analysis  of  ablative 
materials  in  liquid  and  solid 
rocket  motors. 

l 

2. 

Calculation  of  surface  recession  rate 
resulting  from  diffusion  or  kinetically 
controlled  surface  chemical  reactions 
and/or  mechanical  melt  removal. 

Thermal  analysis  of  reentry 
nose  tips. 

3. 

Considers  variable  thermal  properties 
with  a  flexible  physical  model  of 
thermal  conductivity  which  permits  the 
treatment  of  partially  charred  values 
less  than  the  virgin  material  value. 

4. 

Considers  in-depth  and  surface  thermo¬ 
chemical  effects  of  the  pyrolysis  gases 
being  generated  by  both  the  surface  and 
backup  materials. 

ASTHMA 

1. 

Calculation  of  transient  thermal  re¬ 
sponse  of  a  two-dimensional  body. 

Thermal  analysis  of  heat  sink 
materials  in  liquid  and  solid 
rocket  motors  and  reentry  < 

nose  tips. 

2. 

Considers  surface  recession. 

1 

l 

3. 

General  thermochemical  surface  bound¬ 
ary  conditions. 

i 

4. 

Considers  anisotropic  temperature 
dependent  thermal  properties. 

i 

AP.GEIBL 

1. 

Calculation  of  laminar  or  turbulent 
convective  heat  transfer  coefficients 
for  nonablating  wall. 

Transfer  coefficient  analysis 
of  solid  and  liquid  rocket 
motors . 

2. 

Considers  any  real  gas  system  and  axial 
variations  in  wall  temperature  and 
a roe  stream  properties. 

3. 

Solves  the  energy  integral  equation. 

I 

BLIMP 

1. 

Calculation  of  laminar  or  turbulent 
convective  heat  transfer  coefficients 
for  ablating  or  nonablating  wall. 

Transfer  coefficient  analysis 
of  solid  and  liquid  rocket 
motors  and  reentry  nose  cones. 

2. 

Considers  any  real  gas  system  and  axial 
variations  in  wall  temperature,  free 
stream  properties,  and  surface  reces¬ 
sion. 

3. 

Considers  chemical  and  thermal  dlffu-' 
sion  in  the  boundary  layer. 

4. 

Solves  nonsimilar  laminar  and  turbulent 
boundary  layer  equations. 
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TABLE  4-2 


TYPICAL  PROCEDURE  FOR  ANALYZING  ROCKET  MOTOR  NOZZLES 


Sequence 

Purpose 

Applicable 

Computer 

Program 

Result 

Use  of  Results 

1 

Flow  field 

Boundary  layer  edge 
r?ass  flow  rate/unit 
area  as  a  function  of 
distance 

Combined  with  results 
from  Item  2  to  define 
boundary  layer  edge 
properties  as  a  func¬ 
tion  of  distance 

2 

Equilibrium  nozzle 
expansion  analysis 

ACE 

Boundary  layer  edge 
properties  as  a  func¬ 
tion  of  mass  flow 
rate/unit  area 

Combined  with  results 
from  Item  1  as  mentioned 
above .  Information  from 
Items  1  and  2  are  input 
into  ARGEIBL 

3 

Define  Mollier  dia¬ 
gram  and  transport 
properties  for  real 
gas  system 

ACE 

Thermodynamic  and 
transport  property 
maps 

Input  to  ARGEIBL  and 

CMA  or  ASTHMA  computer 
programs 

4 

Define  convective 
heat  transfer 
coefficient 

ARGEIBL 

or 

BLIMP 

Convective  heat  trans¬ 
fer  coefficient  as  a 
function  of  distance 
for  desired  wall  tem¬ 
perature  . 

Input  to  CMA  or  ASTHMA 
computer  program 

5 

Define  radiative 
heat  flux 

Radiative  flux  as  a 
function  of  distance 

Input  to  CMA  or  ASTHMA 
computer  program 

6 

Generation  of  sur¬ 
face  and  in-depth 
thermo chemical 
data 

ACE 

Frozen  edge  enthalpy 
as  a  function  of  tem¬ 
perature  ,  surface  en¬ 
thalpy  as  a  function 
of  surface  temperature 
and  nondimensionalized 
pyrolysis  gas  and  char 
removal  rates  including 
the  effect  of  mechan- 
•ical  melt  removal  if 
applicable,  and  pyrol¬ 
ysis  gas  enthalpy  as  a 
function  of  tempera¬ 
ture 

Input  to  CMA  or  ASTHMA 
computer  program.  Note 
that  the  data  require¬ 
ments  resulting  from  in- 
depth  charring  are  not 
appropriate  to  ASTHMA 

7 

Thermal  response 
predictions 

CMA 

or 

ASTHMA 

In-depth  temperature 
distribution  and  sur¬ 
face  recession  rates 
as  functions  of  time 

Define  insulation  thick¬ 
ness  requirements,  re¬ 
fine  ballistic  predic¬ 
tions,  and  used  as  input 
to  structural  analysis 
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ACE  CODE 


Those  quantities  which  must  be  specified  to  perform  the  thermal  analysis 
using  the  procedures  described  above  are  as  follows: 

•  Propellant  (or  simulator)  elemental  composition 

•  Either  propellant  initial  composition  or  theoretical  flame  temperature 
temperature 

•  Combustion  chamber  operating  conditions  (e.g.,  pressure) 

•  Combustion  product  radiation  properties 

•  Surface  material  elemental  chemical  compositions  and  constituent 
makeup,  if  applicable 

•  Material  virgin  and  char  densities 

•  Material  thermophysical  properties 

Specific  heat 
Thermal  conductivity 
Surface  emissivity 
Heat  of  formation 

Decomposition  kinetic  constants  for 
constituents,  if  applicable 
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Simulator  operating  conditions  and  combustion  product  compositions  have  been 
presented  for  the  five  multiple  pulse  dury  cycles  in  Section  3.  Material  depen¬ 
dent  properties  including  surface  tnermochemical  response  data  are  given  along 
with  the  computed  surface  boundary  conditions  in  Section  4.2. 

4.2  INTERMEDIATE  ANALYSIS  RESULTS  AND  MATERIAL  PROPERTIES 

Analysis  of  material  thermal  response  in  hyperthermal  combustion  product 
environments  requires  the  evaluation  of  convective  and  radiative  boundary  condi¬ 
tion  information  and  the  specification  of  material  thermophysical  properties. 

The  results  of  studies  performed  to  obtain  these  quantities  are  reviewed  in 
this  section. 

The  analysis  of  the  nozzle  flow,  boundary  layer  transport,  and  radiative 
energy  interchange  phenomena  are  covered  in  Section  4.2.1.  These  calculations 
define  the  non-blowing  (i.e.,  non-ablating)  surface  time  dependent  heating 
condition  for  the  five  duty  cycles.  The  ablating  surface  thermochemical  response 
calculations  for  both  the  graphitic  throat  package  and  the  exit  cone  material 
at  A/A*  =  4.2  are  described  in  Section  4.2.2.  Material  thermophysical  properties 
are  given  in  Section  4.2.3. 

4.2.1  Time  Dependent  Non-Blowing  Surface  Boundary  Conditions 

The  time  dependent  non-blowing  surface  boundary  condition  specifications 
may  be  conveniently  divided  into  the  following  categories: 

•  Convective  heating  during  simulator  firings 

•  Radiative  heating  during  simulator  firings 

•  Exposed  surface  boundary  conditions  during  cooldowns 

These  categories  are  discussed  in  Sections  4. 2. 1.1,  4. 2. 1.2,  and  4. 2. 1.3, 
respectively. 

4. 2. 1.1  Convective  Heating  During  Firings 

The  evaluation  of  the  turbulent  convective,  non-blowing  heating  during 
firing  pulses  was  done  utilizing  the  Aerotherm  Real  Gas  Energy  Integral  Boundary 
Layer  (ARGEIBL)  computer  code.  The  uncertainty  in  heat  transfer  coefficient 
using  this  or  any  other  turbulent  boundary  layer  prediction  procedure  has  been 
estimated  to  be  on  the  order  of  +  25  percent  (e.g.,  Reference  7). 

The  ARGEIBL  procedure  requires  the  specification  of  boundary  layer  edge 
properties  as  a  function  of  streamwise  boundary  layer  length.  These  quantities 
were  obtained  from  real  gas,  two-phase  flow  isentropic  nozzle  expansion  solutions 
which  were  computed  using  the  Aerotherm  Chemical  Equilibrium  (ACE)  computer  code. 
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The  expansion  solutions  were  generated  for  the  nominal  simulator  propellant 

composition  given  in  Section  2.  One  solution  corresponded  to  ideal  combustion 

at  700  psia  chamber  pressure  (T  =  6210°R),  and  one  corresponded  to  97  percent 

*  ® 

combustion  efficiency  at  700  psia  (T  =  5830°R) .  The  97  percent  combustion  effi¬ 
ciency  was  simulated  by  doing  equilibrium  solutions  at  a  lower  total  temperature. 

A  comparison  of  the  boundary  layer  edge  static  temperature  ratios  for  the  two 
combustion  efficiencies  is  given  in  Figure  4-3.  As  the  figure  indicates,  the 
effect  of  the  non-ideal  combustion  on  the  expansion  characteristics  is  noticeable. 
The  difference  between  the  two  temperature  ratio  curves  is  due  to  the  greater 
condensed  alumina  formation  at  the  lower  temperature.  The  lower  total  temperature 
also  means  that  the  alumina  solidifies  earlier  in  the  expansion.  Although  Figure 
4-3  indicates  that  the  effect  of  incomplete  simulator  fluid  combustion  should 
be  accounted  for,  errors  in  the  expansion  properties  attributable  to  instantane¬ 
ous  fluctuations  in  combustion  efficiency  or,  equivalently,  in  combustion  product 
composition  are  small  relative  to  the  25  percent  uncertainty  in  the  overall 
boundary  layer  analysis  results. 


ruiUliE  4-3  ISENTROP.C  EXPANSION  TEMPERATURE  RATIOS  FOR  IDEAL  AND 
ACTUAL  PROPELLANT  COMBUSTION  EFFICIENCIES 

*As  used  here,  the  combustion  efficiency  is  defined  as  the  ratio  of  the  actual 
to  the  ideal  characteristic  exhaust  velocity  which  varies  approximately  as 
the  chamber  temperature  to  the  one-half  power. 


The  97  percent  efficiency  solution  was  felt  to  best  characterize  the 
simulator  combustion  product  flow  and  was  used  to  define  boundary  layer  edge 
properties  for  the  ARGEIBL  convective  heating  distribution  computation.  These 
properties  are  shown  as  a  function  of  non-dimensional  mass  flux  in  Figure  4-4. 


FIGURE  4-4  FREE  STREAM  ISENTROPIC  EXPANSION  PROPEP.TIES  * 
COMPUTED  FOR  97  PERCENT  COMBUSTION  EFFICIENCY 


The  enthalpy  values  correspond  to  the  gaseous  component  of  the  two  phase  flow. 

As  discussed  in  Section  4.1,  the  relationship  between  the  expansion  results  and 
the  streamwise  boundary  layer  length  coordinate  is  defined  by  a  nozzle  flow  field 
analysis.  For  the  multiple  pulse  rocket  nozzle  analysis,  the  flow  was  assumed 
to  be  one-dimensional.  Thus  the  ratio  of  local  edge  mass  flux  to  throat  flux 
is  directly  related  to  local  area  ratio  by  mass  conservation,  i.e., 


— —  =  A* /A  =  inverse  of  area  ratio 
pu* 


where  ,.*/A  is  defined  as  a  function  of  streamwise  length  from  the  nozzle  geometry. 
The  radius  and  inverse  area  ratio  for  the  nominal  nozzle  geometry  shown  in 
Figure  2-1  are  shown  as  a  function  of  boundary  layer  length  in  Figure  4-5.  An 
assumption  about  the  flow  over  the  lip  of  the  throat  entrance  cap  was  required 
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FIGURE  4-5  NOMINAL  NOZZLE  GEOMETRY  AS  A  FUNCTION  OF 
STREAMWISE  BOUNDARY  LAYER  LENGTH 


to  define  edge  conditions  in  the  throat  region.  A  separation-recirculation  region 
aft  of  the  throat  entrance  ring  was  expected  but  was  difficult  to  quantify.  For 
simplicity,  the  separation  streamline  shown  in  Figure  4-5  was  assumed.  Based  on 
this  assumption,  streamwise  variations  of  pressure,  velocity,  and  edge  gas  enthalpy 
were  obtained  from  Figures  4-4  and  4-5. 

The  other  information  needed  to  perform  an  ARGEIBL  prediction  are 

1)  Mollier-type  thermodynamic  property  tables  including  Prantl  number  and 

2)  specification  of  the  streamwise  variation  of  gas  phase  enthalpy  at  the  surface 
temperature.  As  explained  in  Section  4.1,  the  Mollier  tables  are  generated  by 
parametric  ACE  code  thermodynamic  state  solutions.  Surface  temperature  variation 
is  estimated  based  on  the  anticipated  material  ablative  response. 

For  the  multiple  pulse  nozzle  analysis,  ARGEIBL  predictions  were  performed 
for  the  anticipated  surface  temperature  (5000°R  isothermal  surface)  and  for  a 
cold  (900°R)  surface  to  assess  the  sensitivity  of  the  heat  transfer  coefficient 
for  a  wide  range  in  surface  temperature.  The  two  transfer  coefficient  distribu¬ 
tions  are  shown  in  Figure  4-6.  Based  upon  the  nozzle  throat  values  of  the  con¬ 
vective  heat  transfer  coefficient,  the  two  solutions  indicate  that  the  coefficient 
varies  as  temperature  to  the  -0.22  power.  This  means  that  a  20  percent  error 
in  estimating  surface  temperature  (1000°R)  amounts  to  only  a  4  percent  error  in 
heat  transfer  coefficient. 
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FIGURE  4-6  ARGF.IRL  HEAT  TRANSFER  COEFFICIENT  DISTRIBUTION  FOR  MULTIFLE  PUL'SE 
NOZZLE  COI. FI GL* PAT  ION.  AN8  3066  SIMULATOR  PROPELLANT 


To  this  point,  the  convective  heating  evaluation  was  performed  for  the 
nominal  rocket  simulator  operating  conditions  (propellant  composition  from 
Table  2-3  and  a  combustion  chamber  pressure  of  700  psia) .  Turbulent  heating 
considerations  (Reference  8)  indicate  that  in  a  rocket  nozzle,  the  heat  transfer 
coefficient  varies  as  the  0.8  power  of  chamber  pressure.  To  account  for  pressure 
changes  during  the  various  firing  pulses,  the  following  relationship  was  used: 


peueCH  "  (K) 
o 


I  foT 

\  700/  ( 


peueCH 


ARGEIBL 


(4-1) 


where 


PGueCH  -  non-blowing  convective  heat  transfer  coefficient  appli- 

°  cable  to  a  recovery  enthalpy  driving  potential 

PQ  -  instantaneous  measured  chamber  pressure  in  psia  units 

(see  Section  2  and  Appendix  A) 


peueCI! 


ARGEIBL 


predicted  ARGEIBL  heat  transfer  coefficient  at  the 
nozzle  location  of  interest 


K 


-  correlation  factor  (evaluated  in  Reference  9  as  0.75) 


The  variation  of  heat  transfer  coefficient  with  time,  plus  the  applicable  local 
recovery  enthalpy,  defines  the  convective,  time  dependent  boundary  condition  dur¬ 
ing  the  firing  pulses. 
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In  applying  the  convective  heating  boundary  condition  as  a  function  of 
time  to  the  PG  coated  nozzle  throat,  consideration  had  to  be  given  to  the  protec¬ 
tion  of  the  upstream  section  of  the  coating  by  the  overlapping  portion  of  the 
MXC  313  carbon  phenolic.  This  protection  of  the  PG  coating  was  taken  into  account 
by  assuming  zero  convective  heating  to  that  portion  of  the  coating  which  was 
covered  by  the  MXC  313  carbon  phenolic  until  such  time  as  the  phenolic  carbon 
had  recessed  and  exposed  the  coating.  In  addition,  convective  heating  in  the 
separation  region,  which  is  shown  in  Figure  4-5  for  the  initial  configuration, 
was  assumed  to  be  negligible.  As  the  MXC  313  carbon  phenolic  recedes,  the  axial 
length  of  this  separation  region  decreases. 

4. 2. 1.2  Radiation  Heating  During  Simulator  Firings 

In  both  the  CMA  and  the  ASTHMA  codes,  the  surface  energy  balance  takes 
the  following  form: 

^conv  +  ^chem  prod  +  ^rad  in  ~  ^rad  out  ^cond  ®  (42) 

These  flux  terms  account  for  the  energy  transport  phenomena  shown  in  Figure  4-1. 
This  section  is  concerned  with  the  two  radiation  terms.  The  radiant  flux  absorbed 
.by  the  surface  (4ra<j  ^n)  is  defined  as  the  incident  radiation  multiplied  by  the 
absorptance  of  the  wall.  The  radiant  flux  emitted  by  the  surface  ^  out)  is 
defined  as  that  emitted  by  a  black  body  multiplied  by  an  emittance  factor  and 
by  a  view  factor.  Therefore  the  net  radiation  to  the  surface  (<5net  ra<j)  is 
defined  by 


L  .  =  ad.  -  Fc  foT* 

net  rad  w’mc  rad  w  w 


(4-3) 


where 

o  , e  -  absorptivity  and  emissivity,  respectively,  of  wall  material 
(<*w  and  ew  are  equal  by  Kirchoff's  Law  for  a  gray  body) 

a  -  Stefan-Boltzmann  constant 
T^  -  wall  temperature 
F  -  view  factor 

In  order  to  model  the  radiation  transport  between  a  particle  laden  stream 
and  the  nozzle  surface,  it  was  assumed  that  the  stream  was  optically  thick  and 
that  the  particles  and  wall  exchange  radiant  energy  as  if  they  were  two  parallel 
plates.  In  this  way  multiple  reflections  between  the  wall  and  stream  were  taken 
into  account.  In  addition,  it  was  assumed  that  both  the  stream  and  wall  behave 


as  gray  bodies  and  that  they  emit  and  reflect  radiant  energy  diffusely.  Based 
on  the  above  assumption,  the  net  radiant  heat  flux  relation  is  given  in  as: 


^net  rad  ~  eeff*oTs 


oO 


w 


where 


eeff  (e3^ective  emissivity) 


es  -  particle  laden  stream  emissivity 
Tg  -  local  stream  static  temperature 


(4-4) 


(4-5) 


To  define  the  net  radiation  heat  flux  using  Equation  (4-4)  the  local 
static  stream  temperature  and  the  effective  emissivity  had  to  be  evaluated. 

The  stream  temperature  was  obtained  from  the  isentropic  expansion  data  presented 
previously  in  Figure  4-3.  At  the  nozzle  throat  and  upstream  thermocoupied  plug 
analysis  prediction  locations,  the  stream  temperatures  are  5000  and  3710°F, 
respectively.  The  effective  emissivity  was  evaluated  from  Equation  (4-5)  using 
a  value  of  0.9  for  the  wall  emissivity  and  using  the  following  relationship 
(Reference  8)  for  the  emissivity  of  the  particle  laden  combustion  products. 


es  =  1  -  exp 


(4-6) 


where 

C  -  experimentally  determined  constant  (0.808  for  the  propellant 
being  considered  m  this  analysis) 

n  -  percentage  of  aluminum  loading  {=  15.5  for  ANB  3066  simulation 
propellant) 

p  -  local  density  of  propellant  combustion  species  (lb/ft3) 

D  -  local  beam  length,  usually  taken  as  the  diameter  (in.) 

At  the  nozzle  throat  and  exit  cone  analysis  prediction  locations,  the 
effective  emissivities  calculated  from  Equation  (4-5)  were  0.30  and  0.09,  respec¬ 
tively,  for  a  nominal  chamber  pressure  of  700  psia.  However,  a  distribution  of 
the  effective  emissivity  with  streamwise  length  was  required  for  the  two-dimen¬ 
sional  analysis  of  the  nozzle  throat  insert  for  nozzle  no.  8.  This  distribution 
is  shown  in  Figure  4-7  for  each  of  the  four  pulses  which  comprised  the  multiple 
pulse  duty  cycle  for  this  nozzle.  This  figure  shows  the  increase  in  effective 
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FIGURE  4-7,  VARIATION  OF  EFFECTIVE  EMISSIV ITY  WITH  STREAMHISE  LFNGTH 

emissivity  with  the  recession  of  the  MXC  313  carbon  phenolic  entrance  cap. 

Average  values  of  local  pressure  (or  area  ratio)  during  a  particular  pulse  were 
used  for  predicting  the  effective  emissivity  as  the  accuracy  of  the  stream 
emissivity  correlation  (Equation  (4-6))  did  not  warrant  a  detailed  analysis. 

In  conclusion,  the  overall  effect  of  radiation  as  compared  to  convective 
heat  transfer  during  the  firing  pulses  is  small.  For  example,  during  the  first 
pulse  of  the  nozzle  no.  8  firing,  the  net  radiation  at  an  exit  cone  area  ratio 
of  4.2  amounted  to  only  4  percent  of  the  total  energy  transferred  to  the  nozzle 
surface. 

4. 2. 1.3  Exposed  Surface  Boundary  Condition  During  Cooldown 

At  motor  burnout,  convective  heating  of  the  nozzle  surface  is  eliminated 
which  simplifies  the  surface  energy  balance  given  by  Equation  (4-2)  to  the 
following  form: 


^rad  in  ^rad  out  ^cond  -  ® 


(4-7) 
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In  applying  this  equation  during  cooldown ,  it  was  assumed  that  the  nozzle 
surface  was  an  isothermal  gray  body.  This  assumption  implies  there  is  no  net 
radiant  heat  transfer  between  the  surface  under  consideration  and  the  rest  of 
the  nozzle  interior.  The  only  net  transfer  of  radiant  energy  occurs  between  the 
section  of  wall  being  analyzed  and  the  environment  beyond  the  exit  plane  of  the 
nozzle  exit  cone.  This  flux  is  given  by 

<*net  rad  =  Few°  <Te  “  <4‘8) 

where 

Tg  *  environment  temperature  exterior  to  the  exit  cone 
Neglecting  T*  in  comparison  to  T^,  this  equation  becomes 

W  rad  *  -  FEw0Ti  <4-9> 

Equation  (4-9)  was  used  to  model  the  radiative  heat  flux  during  cooldown. 

In  evaluating  <5net  ra^  at  the  nozzle  throat  and  exit  cone  analysis  loca¬ 
tions  from  Equation  (4-9)  a  value  of  0.9  was  used  for  the  wall  emissivity  and 
values  of  0.011  and  0.310  were  used  for  the  view  factor  at  the  nozzle  throat  and 
exit  cone  locations  respectively.  The  view  factors  were  calculated  using 
the  procedures  presented  in  References  9  and  10. 

A  brief  analysis  was  performed  to  study  the  sensitivity  of  the  predicted 
cooldown  temperatures  on  the  magnitude  of  the  view  factor  for  the  upstream 
thermocoupled  plug.  For  this  sensitivity  analysis,  one  prediction  was  made  using 
the  view  factor  mentioned  above  (0.310)  and  one  was  made  assuming  a  view  factor 
of  zero.  The  temperatures  were  predidW.  for  the  cooldown  following  pulse  no.  1 
for  nozzle  no.  8,  and  these  temperatures  are  shown  in  Figure  4-8.  Based  on  these 
results  two  conclusions  are  apparent,  namely 

•  the  predicted  temperatures  during  cooldown  are  a  weak  function  of  the 
radiation  view  factor 

•  the  predicted  temperatures  during  cooldown  at  the  TC  no.  1  location 
for  either  value  of  the  radiation  view  factors  are  considerately  below 
the  measured  temperatures 

Even  though  the  finite  value  for  the  radiation  view  factor  causes  the  cooldown 
temperatures  to  be  under  predicted,  it  best  represents  the  radiation  phenomena 
during  cooldown,  and,  thus,  it  was  used  for  all  of  the  multiple  pulse  nozzle 
analyses.  The  reason  for  the  temperatures  being  under  predicted  during  cooldown 
is  presented  in  Section  4.3. 


During  the  firing  pulses,  the  nozzles  were  exposed  to  a  high  temperature 
chemically  reacting  environment  which  caused  chemical  corrosion  at  the  surface 
of  all  of  the  components  and  caused  in-depth  material  degradation  and  off-gassing 
for  the  charring  ablating  components.  In  order  to  account  for  the  complicated 
boundary  layer  thermochemical  phenomena  occuring  at  the  nozzle  surface  in  such  an 
environment,  the  ACE  code  described  in  Section  4.1  was  utilized  to  generate 
parametric  sets  of  surface  thermochemical  solutions  suitable  for  input  to  the 
CMA  or  ASTHMA  thermal  response  analysis  procedures.  The  surface  thermochemistry 
results  computed  for  the  throat  package  are  presented  in  Section  4. 2. 2.1  and 
those  for  the  exit  cone  location  are  presented  in  Section  4. 2. 2. 2. 
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4. 2. 2.1  Graphite  Surface  Thermochemistry  Solutions 

The  surface  thermochemical  response  of  graphitic  materials  such  as  ATJ 
and  pyrolytic  graphite  (PG)  is  controlled  by  the  finite  (as  opposed  to  equilibrium) 
reaction  rates  of  the  surface  carbon  atoms  with  the  HjO,  C02  and  H2  components 
of  the  solid  propellant  simulator  combustion  products.  The  kinetically  controlled 
reaction  rate  models  developed  under  Air  Force  Contract  F04611-69-C-0081  were 
utilized  to  compute  the  surface  response  of  the  graphite  components  in  the  multi¬ 
ple  pulse  nozzle  throat  package.  The  solutions  were  performed  using  the  ACE  code. 

The  kinetic  response  model  incorporates  a  Langmuir  type  surface  reaction 
inhibition  expression  into  an  Arrenhius  temperature  dependent  rate  relation. 

The  rate  expression  and  reaction  kinetic  rate  constants  are  given  in  Table  4-3. 


TABLE  4-3 

KINETICALLY  CONTROLLED  REACTION  RATE  EXPRESSIONS  USED  FOR  PG  LAYER 
AND  ATJ  GRAPHITE  SURFACE  THERMOCHEMISTRY  SOLUTIONS11 


Material 

Pre-Exponent  Factors 
(units  to  give  below) 

Activatibn  Energies 
(ca 1/mole) 

D1 

°2 

E1 

E2 

Layer  PG 

12.5 

0.77 

40,000 

55,500 

Edge  PGb 

4.12  x  10' 

1.03  X  10‘ 

65,500 

129,500 

ATJ  (Bulk)  Graphite 

1.98  X  10S 

4.94  x  107 

i 

1 

*„  _  B1cxp(-E1/RTW) 

xc*  “  u.  ♦  ph  0  *  P|l  +  pc0  *  pco>  rrr 


7  .  !C0P"A  (  .  PC0  s 
>HCt>  y  H2°  KHj0  J  yPC02  KC0^ 


fl2CRp(-K,/RTw) 

T17"*  PlIjO  *  pco2  +  PC0T 


(  Pc?*2 
1P!,2  ‘  KH 


,  rcoics  of  rcaction/f L:scc 


wo  edge  PG  solutions  were  performed;  constants  included  for  completeness. 


The  partial  pressure  quantities  were  evaluated  in  the  ACE  code  using  the 
film  coefficient,  multicomponent  unequal  diffusion  boundary  layer  model.  Gaseous 
species  adjacent  to  the  graphite  surface  were  allowed  to  be  in  chemical  equil¬ 
ibrium  with  each  other  but  not  with  the  surface.  In  order  to  allow  diffusion 
of  only  ♦•he  gaseous  components  of  the  frecstream  fluids,  an  effective  boundary 
layer  eage  gas  composition  was  defined  by  removing  those  elemental  amounts  tied 
to  the  condensed  alumina  from  the  ANB  3066  composition.  The  relative  elemental 
composition  with  and  without  alumina  are  given  in  Table  4-4 . 
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"‘ABLE  4-4 

CHEMICAL  ELEMENTAL  COMPOSITIONS  POP  THROAT  AND  EXIT  CONE 
SURFACE  THERMOCHEMISTRY  SOLUTIONS* 


Element 

Boundary  Layer  Edge 

PE  i  ATJ 
Surface 

FM5055A  Carbon 

Phenolic 

■\ 

I 

NO 

a12°3*w 

Allowed" 

Pyrolysis 

Gas 

Char 

Hydrogen 

3.460 

3.460 

— 

10.597 

— 

Carbon 

0.905 

0.905 

8.326 

5.083 

8.326 

Nitrogen 

0.624 

0.624 

— 

... 

— 

Oxygen 

2.3B2 

1.534 

— 

1.767 

Aluminum 

0.583 

0.018 

— 

— 

- - 

Chlorine 

0.651 

0.551 

— 

— 

— 

aAl)  quantities  in  gr  atomn/100  gr  of  mixture 

bhass  of  AljOj*  removed  per  muss  of  gas  phase  remaining  =  0.404 


Surface  maps  were  generated  for  two  locations  for  both  the  ATJ  throat 
retainer  ring  and  the  layer  oriented  pyrolytic  graphite  throat  insert.  One 
of  the  pyrolytic  graphite  maps  was  generated  for  the  throat  plane  conditions  and 
used  for  all  one-dimensional  throat  response  calculations.  The  conditions  for 
the  other  three  maps  were  defined  so  that  streamwise  variations  in  surface  ther¬ 
mochemical  response  could  be  accounted  for  in  the  two-dimensional  throat  package 
response  calculation  (nozzle  no.  8  only).  The  area  ratio,  local  static  pressure 
and  temperature,  local  edge  gas  recovery  enthalpy,  and  the  local  mass  transfer 
coefficient  values  for  which  the  kinetically  controlled  response  calculations 
were  performed  are  given  in  Table  4-5.  The  non-dimensional  ablation  rate  maps 

TABLE  4-5 

SUMMARY  OF  THROAT  AND  EXIT  CO.’E  CONDITIONS  FOR  WHICH 
SURFACE  THEP.MOCKEMI  STRY  SOU/ 'IONS  WERE  OBTAINED 


■ 

Edge  Static 

Conditions 

Nominal  Transfer  Coefficient 

Recovery . 
Enthalpy0 

(Btu/lbm) 

Pressure 

(atm) 

Temperature 

(°ro 

Meat*3 

(Ib/ft'soc) 

Mass0 

(Ib/Ct Jsec) 

-1.7 

PG 

44.0 

3208. 

" 

0.67 

0.53 

*.P7. 

Throat 

PG 

27.6 

3029. 

0.82 

0.65 

361. 

1.4 

ATJ 

13.5 

2756. 

0.53 

0.42 

324 . 

2.0 

ATJ 

f  .5 

2321. 

0.36 

23 

26.’. 

4.2 

FM5055A 

2.4 

2313. 

0.17 

0.1- 

_ ill: _ 

NegAtive  nrcbor  indicates  subronlr  c rc.'.  ratio 
bHoat  transfer  coefficient  equals  0.71  tires  AROEIRL  result 
^ratio  o'  mars  to  heat  coefficient  equals-  C.734 
^Recovery  enthalpy  equals  h3TfVrlc  *  ?r'^* (u’/7) :  T  »  0.504 
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FIGURE  4-5  KINETICALLY  CONTROLLED  RESPONSE  OF  PYROLYTIC  AND  ATu 
GRAPHITES  TO  AN3  306G  PROPELLANT  SIMULATION 

for  each  location  are  shown  in  Figure  4-9.  Also  shown  in  Figure  4-9  are  the 
peak  surface  temperature  points  attained  during  the  ASTHMA  throat  package 


thermal  analysis. 


4  •  2 .  c .  « 


ctrbon  Cloth  Phenolic  Exit  Cone  Thermochemistry  Solution 


Diffusion  controlled  surface  thermochemistry  state  solutions  were  com¬ 
puted  using  the  ACE  computer  coue  for  conditions  corresponding  to  an  area  ratio 
of  4.2  in  the  exit  cone  of  the  multiple  pulse  nozzles.  These  solutions  para¬ 
metrically  defined  the  surface  response  analysis.  The  surface  thermochemistry 
solutions  utilized  1)  the  same  multicomponent,  unequal  diffusion  coefficient 
boundary  layer  model  that  was  used  for  the  kinetically  controlled  throat  response 
analyses  (Section  4. 2. 2.1),  2)  the  effective  edge  gas  composition  given  in 

Table  4-3  as  the  change  in  the  freestream  alumina  between  the  throat  an!  the 
A/A*  *  4.2  locations  was  negligible,  and  3)  the  boundary  layer  edge  conditions 
given  in  Table  4-4.  Furthermore,  since  the  carbon  cloth  phenolic  depolymer izes 
in  depth  at  varying  rates,  the  effect  of  different  relative  amounts  of  pyrolysis 
gas  on  the  surface  thermochemical  response  was  evaluated.  The  computed  non- 
dimensional  ablation  rates  for  parametrically  changing  values  of  non-dimensional 
pyrolysis  off  gas  rate  are  shown  in  Figure  4-10. 


FIGURE  4-10  DIFFUSION  CONTROLLED  SURFACE  RESPONSE  OF  FK5055A  CARSON  PHENOLIC  TO 

AN33065  PROPELLANT  SIMULATION  FCR  VARIOUS  NON-DIMENSIONAL  PVRClVSIS  GAS  SATES. 

4.2.3  Thermal  Properties 

A  brief  literature  survey  was  conducted  to  define  the  most  applicable 

thermal  properties  to  be  used  in  predicting  the  thermal  performance  of  the 

FM  5055A  carbon  cloth  phenolic.  This  survey  revealed  that  the  desired  data  could 

be  obtained  from  either  Reference  2  or  Reference  11.  The  major  descrepency 

between  the  data  presented  in  these  references  is  for  the  virgin  and  char  thermal 

conductivities.  In  Reference  2,  virgin  and  char  thermal  conductivities  are 

* 

presented  for  MX  4926  while  in  Reference  11,  a  model  is  presented  for  calcula¬ 
ting  the  virgin  and  char  thermal  conductivities  for  either  MX  4926  or  FM  50f>5A. 

The  thermal  conductivity  results  from  these  references  are  shown  in  Figure  4-11 
with  the  properties  labeled  MX  4926  obtained  from  Reference  2  and  with  the  proper¬ 
ties  labeled  FM  5055A  obtained  from  Reference  11.  Since  virgin  material  thermal 
conductivity  above  20008R  has  very  little  significance  because  the  material  is 
generally  approaching  the  fully  charred  state  at  this  temperature,  the  virqin 
and  char  thermal  conductivities  from  the  two  references  for  a  90  degree  layup 
angle  are  in  close  agreement.  In  addition  the  virgin  material  properties  from 


* 

MX  4926  and  FM  hOSSA  are  similar  materials  which  arc  fabricated  by  different 
manufacturers. 


the  two  sources  are  in  close  agreement  for  the  0  degree  layup  angle  material. 
However,  Figure  4-lla  shows  that  a  major  discrepency  exists  for  the  char  thermal 
conductivity  from  the  two  references  for  the  0  degree  layup  material.  Since 
the  layup  angle  for  the  FM  5055A  exit  cone  under  consideration  is  15  degrees, 


CHAR 

1 

v  t  f)P  t  H 

V inuiH 

/ 

A 

TEMPERATURE ,  DEGREES  R 
0  DEGREE  LAYUP 


VIRGIN 


■FM  5055A 


MX  4925 


TEMPERATURE,  DEGREES  R 
b)  90  DEGREE  LAYUP 

FIGURE  4-1!  THERMAL  CONDUC  i  IV I TY  OF  FM  5055A  AND  <X  4926 
CARBON  CLOTH  PHENOL  ICS 

the  discrepency  in  the  char  thermal  conductivities  at  a  0  degree  layup  angle  had 
to  be  resolved.  As  a  first  step,  CMA  analyses  we re  performed  for  pulse  no.  1 
of  nozzle  no.  4  using  both  sets  of  thermal  conductivity  properties,  and  the  pre¬ 
dicted  temperatures  from  the  analyses  are  compared  with  the  measured  in-depth 
temperatures  in  Figure  4-12.  In  performing  these  analyses,  the  input  information 
{Section  4.1)  other  than  the  virgin  and  char  thermal  conductivities  was  obtained 


FIGURE  4-12  COMPARISON  OF  MEASURE!)  AND  PREDICTED  TEMPERATURES  FOR  DIFFERENT  CARBON 
CLOTH  PHENOLIC  THERMAL  PROPERTIES 


from  Reference  2  and  a  portion  of  this  information  is  summarized  in  Table  4-6. 

TABLK  4-6 


PHYSICAL  AND  CHEMICAL  PROPERTIES 
TOR  "K  jOfO.1  CAl’WfJ  CLOTH  phenolic 


Nominal  Density,  Ib/Ct ’ 

89.4 

Ilorrinal  n  ’/action 

0.345 

Renin  Residual 

0.40 

Virgin  Material  Heat  of  Karma lion,  utu/lb 

-379 

Assumed  Renin  Elemental  Formula 

C61!C° 

Reinforcement  Elemental  Formula 

C 

1 _ 

In  addition,  the  effect  of  material  layup  angle  on  the  thermal  conductivity  was 
accounted  for  by  using  the  following  relationship: 


(4-9) 


where 


k0°'  k90e 


--  0 


thermal  conductivity  for  the  0°  and  90°  layup  angles, 
respectively 

layup  angle 
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Careful  examination  of  Figure  4-12  shows  that  the  predicted  temperatures 
using  the  MX  4926  properties  are  in  good  agreement  with  the  measured  temperatures 
for  temperatures  above  1200°R  but  that  the  predicted  temperatures  using  FM  5055A 
properties  provide  better  agreement  for  temperatures  below  1200°R.  A  review  of 
the  comparisons  between  the  measured  and  predicted  temperatures  presented  in 
Reference  2,  which  were  for  MX  4926  thermal  properties,  also  showed  that  the 
material  temperatures  were  overpredicted  for  temperatures  less  than  12008R, 

Based  on  the  temperature  comparisons  presented  in  Figure  4-12  and  Reference  2, 
a  set  of  modified  FM  5055A  thermal  conductivities  were  defined  for  a  0  degree 
layup  angle.  These  properties  which  combine  the  best  features  from  the  MX  4926 
cind  FM  5055A  properties  are  presented  in  Figure  4-13. 


The  predicted  temperatures  using  these  modified  properties  are  also  presented  in 
Figure  4-12  and,  as  expected,  provide  better  overall  agreement  with  the  measured 
temperatures.  As  a  further  check  on  the  modified  FM  5055A  thermal  properties, 
option  2  CMA  analyses  were  performed  using  TC  no.  2  of  pulses  no.  1  and  2 
for  nozzle  no.  8  as  an  impressed  boundary  condition.  For  the  MX  4926,  FM  5055A, 
and  modified  FM  5055A  thermal  properties,  the  predicted  in-depth  temperatures 
at  the  TC  no.  3  location  are  compared  with  the  measured  temperatures  in  Figure  4-14. 


FIGURE  4-14  COMPARISON  OF  PREDICTED  AMD  MEASURED  1N-0EPTH  TEMPERATURE  USING 
OPTION  Z  CMA 

The  temperatures  predicted  using  the  modified  FM  5055A  properties  are  in  better 
agreement  with  the  measured  data  than  those  using  the  MX  4926  properties.  How¬ 
ever,  the  slope  of  the  measured  temperatures  for  pulse  no.  2  and  the  subsequent 
cooldown  was  predicted  better  by  the  temperatures  using  the  MX  4926  thermal 
conductivities. 

Based  on  the  above  analyses,  the  modified  FM  5055A  thermal  conductivities 
were  considered  the  most  realistic  and  were  used  in  predicting  the  thermal  perfor¬ 
mance  of  the  FM  5055A  carbon  cloth  phenolic  at  a  supersonic  area  ratio  of  4.2. 

The  thermal  properties  used  for  the  PG  coating,  the  ATJ  and  AGSR  graphites,  the 
GA  carbon,  and  the  RPD  150  asbestos  phenolic  materials  in  the  nozzle  throat 
analysis  are  presented  in  Figure  4-5  and  Tables  4-7  and  4-8, 

4-r  i 
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TABLE  4-7 


SPECIFIC  HEAT  AMO  DENSITIES  FOR  PYROLVTIC  GRAPHITE, 
ATJ  GRAPHITE,  AGSR  GRAPHITE,  AND  GA  CARBON 


ATJ  (p 
and 
to  » 

*■  108  lb/ft’) 

GA  Carbon 

96.6  lb/fts) 

Pyrolytic  Graphite 
(p  =  137.3  lb/ft1) 

AGSR  Graphite 
(p  -  96.1  lb/ft’) 

Temp 

Specific  Heat 

Temp 

S-ccific  Heat 

Temp 

Specific  Heat 

<°R) 

(Btu/lb’F) 

(°R) 

(Btu/lb°F) 

( °F) 

(Btu/lb'F) 

492 

.170 

525 

.240 

.220 

720 

.250 

760 

.276 

.323 

1080 

.336 

1060 

.313 

.390 

1900 

.415 

1650 

.409 

KJjTjljW 

.48) 

2160 

.440 

2190 

.445 

4000 

.504 

2520 

.460 

2650 

.480 

5000 

.520 

2880 

.475 

3060 

.498 

6000 

.530 

3240 

.483 

4460 

.520 

3960 

.505 

7460 

.526 

4680 

.514 

5400 

i _ 

.524 

TADLi:  4-8 

THERMAL  PROPERTIES  USED  FOR  RPD  150  ASBESTOS  PHENOLIC* 
(p  =  90.6  lb/ft1) 


Temp 

CR) 

Specific  He.it 
(Htu/lbT) 

Thermal 
Conductivity 
(Btu/f t  scc'R) 

492 

.170 

.825  x  10"’ 

720 

.250 

.825  x  10"k 

1000 

.336 

.825  x  10"k 

1440 

.385 

.850  i:  10“ 

1800 

.415 

.850  x  10“ 

2160 

.440 

.870  x  10"fc 

2520 

_ 

.460 

.900  x  10"“ 

_ _  _ 

Since  the  ASTHMA  Program  docs  not  account  fer 
in-depth  charring,  the  virgin  material  properties 
have  been  modified  to  take  this  into  account. 


4.3  COMPARISON  OF  COMPUTED  AND  MEASURED  MATERIAL  PERFORMANCE  DATA 

Based  on  the  measured  multiple  pulse  duty  cycle  operating  conditions, 
one  dimensional  transient  ablation  and  thermal  response  analyses  were  performed 
for  both  the  upstream  thermocoupled  plug  location  (exit  cone  A/A*  =  4.2)  and  the 
throat  plane  of  each  of  the  five  nozzles.  For  nozzle  no.  8,  a  two-dimensional 
analysis  was  done  for  the  graphite  throat  package  and  its  asbestos  phenolic  backup 
material.  Comparison  of  the  results  from  the  one-dimensional  solution  to  those 
from  the  two-dimensional  solution  for  this  nozzle  gave  an  accurate  definition  of 
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the  differences  between  these  two  types  of  solutions.  These  differences  at  the 
back  surface  of  the  AGSR  substrate  were  applied  to  the  one-dimensional  solutions 
for  the  other  four  nozzles  in  order  to  arrive  at  approximately  two-dimensional 
solutions.  The  thermal  performance  predictions  for  the  nozzle  throat  package 
and  for  the  FM  5055A  nozzle  exit  cone  are  compared  to  the  measured  material  thermal 
performance  in  Sections  4.3.1  and  4.3.2,  respectively. 

4.3.1  Throat  Package  Response  Analysis  Results 

A  two-dimensional  thermal  response  analysis  including  the  effects  of  sur¬ 
face  ablation  was  performed  for  the  throat  package  of  nozzle  no.  8.  This  solution 
provided  an  accurate  evaluation  of  the  temperature  field  within  the  throat 
package  for  the  entire  four  pulse  duty  cycle.  The  nodal  layout  for  this  ASTHMA 
2-D  analysis  as  well  as  the  locations  for  the  two  thermocouples  is  shown  in 
Figure  4-16.  Part  of  the  asbestos  phenolic  backup  material  is  not  included  in 


FIGURE  4-16  THROAT  PACKAGE  NODAL  LAYC  T  fOR  2-D  ANALYSIS 
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the  nodal  grid  as  it  does  not  influence  the  thermal  response  of  the  throat  pack¬ 
age.  The  surface  thermochemistry  maps  shown  previously  in  Figure  4-9  were 
defined  for  surface  locations  corresponding  to  grid  system  axial  distances  (see 
Figure  4-16)  of  0.5,  the  throat,  1.5  and  2.5  inches. 

To  make  maximum  utilization  of  the  2-D  analysis  for  nozzle  no.  8,  a  1-D 
analysis  was  performed  at  the  nozzle  throat  plane  using  the  CMA  code.  By  com¬ 
paring  the  surface  recession  and  the  surface  and  in-depth  transient  temperatures 
from  both  solutions,  the  adequacy  of  the  1-D  analysis  can  be  defined.  Since  the 
surface  energy  balance  formulation  in  the  CMA  and  ASTHMA  codes  are  identical,  the 
same  throat  surface  response  map  could  be  utilized  in  each  solution.  The  convec¬ 
tive  and  radiative  surface  boundary  conditions  could  also  be  specified  identically 
so  that  the  only  difference  between  the  CMA  and  ASTHMA  solutions  at  the  throat 
plane  is  the  2-D  heat  conduction  effects.  Since  the  materials,  geometries,  and 
overall  heat  flux  levels  are  similar  for  all  five  nozzle  firings,  the  differences 
in  the  temperatures  predicted  by  the  two  types  of  solutions  for  nozzle  no.  8  can 
be  applied  to  the  temperatures  predicted  by  the  1-D  throat  analyses  for  nozzles 
no.  4,  5,  6,  and  7  in  order  to  obtain  approximate  2-D  predicted  temperatures.  By 
comparing  these  temperatures  to  the  measured  temperatures  at  the  back  wall  of  the 
throat  insert,  the  relative  accuracies  of  the  two  predicton  procedures  can  be 
evaluated  for  each  of  the  five  duty  cycles.  Comparisons  of  computed  one-  and 
two-dimensional  thermal  response  results  are  compared  to  the  measured  back  wall 
temperature  histories  in  Section  4. 3. 1.1.  The  throat  plane  surface  recession 
calculations  are  compared  to  data  in  Section  4. 3. 1.2. 

4.3.1..  In-Depth  Temperature  Comparisons 

The  predicted  transient  temperatures  using  the  ASTHMA  2-D  analysis  are 
compared  in  Figure  4-17  with  the  temperatures  measured  at  the  rear  surface  of 
the  AGSR  for  nozzle  no.  8.  As  shown  by  this  figure,  the  predicted  temperatures 
which  take  two-dimensional  effects  into  account  are  in  excellent  agreement  with 
the  measured  data  as  the  difference  between  the  measured  and  predicted  2-D  temp¬ 
eratures  is  less  than  100°F  for  the  entire  nozzle  no.  8  duty  cycle  (4  pulses). 
Figure  4-18  shows  the  measured  and  predicted  temperature  histories  for  the  thermo¬ 
couple  at  the  back  wall  of  the  ATJ  throat  retainer  ring.  Again,  the  temperatures 
predicted  by  the  two  dimensional  analysis  are  in  excellent  agreement  with  the 
measured  temperature  as  the  maximum  difference  throughout  all  four  pulses  is  less 
than  130°F . 

The  close  agreement  between  the  2-D  predicted  and  the  measured  temperature 
histories  indicates  that  the  overall  energy  content  of  the  throat  package  has  beer, 
correctly  computed.  Even  though  this  energy  content  is  the  integrated  effect 
of  the  surface  convective  and  radiative  heat  fluxes  during  the  entire  duty 
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cycle,  the  agreement  between  the  predicted  and  measured  temperatures  does  sub¬ 
stantiate  the  analysis  procedures  used  to  define  the  environmental  boundary  con¬ 
ditions  at  the  nozzle  surface. 

The  results  of  the  one-dimensional  CMA  analysis  for  nozzle  no.  8  are  also 
shown  in  Figure  4-17.  The  one-dimensional  predictions  closely  parallel  the  two- 
dimensional  ASTHMA  predictions  but  are  consistently  low  by  as  much  as  550°F. 

This  result  was  expected  because: 

•  The  path  of  thermal  energy  transport  in  the  1-D  analysis  is  radial 
and,  thus,  all  of  the  heat  must  be  conducted  through  the  PG  coating 
which  has  a  thermal  conductivity  about  15  times  smaller  than  for  both 
the  AGSR  substrate  and  the  ATJ  throat  retainer  ring  materials. 

•  The  path  of  thermal  energy  transport  in  the  2-D  analysis  is  both  radial 
and  axial  and,  thus,  because  of  the  lower  thermal  resistance  of  the 

ATJ  throat  retainer  ring,  heat  enters  this  ring  from  the  nozzle  environ¬ 
ment  and  is  conducted  axially  into  the  AGSR  substrate;  this  axial  flow 
of  heat  causes  the  temperatures  predicted  by  the  2-D  analysis  to  be 
higher  than  those  predicted  by  the  1-D  analysis. 

The  remaining  four  nozzles  (nozzles  no.  4,  5,  6,  and  7)  were  analyzed 
using  only  the  one-dimensional  CMA  analysis.  To  estimate  the  expected  two- 
dimensional  analysis  results  for  these  nozzles,  a  simple  model  was  formulated 
to  account  for  the  axial  heat  conduction.  In  this  analysis,  the  amount  of  thermal 
energy  which  arrived  at  the  throat  thermocouple  was  assumed  to  come  from  two 
separate  fluxes.  One  flux  was  that  which  penetrates  the  PG  layer  and  the  second 
flux  is  that  which  is  conducted  axially  through  the  ATJ  graphite.  This  second 
flux  was  computed  in  the  two-dimensional  analysis  but  not  in  the  one-dimensional 
analysis.  These  heat  fluxes  were  correlated  with  the  temperature  rise  of  the 
thermocouple  at  the  backside  of  the  AGSR  at  the  nozzle  throat  plane.  The 
temperature  rise  due  to  heat  flux  through  the  PG  coated  throat  is  predicted  by 
the  one-dimensional  CMA  analysis.  The  additional  temperature  rise  resulting 
from  the  two-dimensional  temperature  prediction  is  a  result  of  axial  heat  conduc¬ 
tion.  The  temperature  difference  between  the  two  prediction  techniques  at 
the  backside  of  the  AGSR  are  given  in  Figure  4-19  as  a  function  of  the  local 
temperature  predicted  by  the  one-dimensional  analysis.  Using  the  correlation  pre¬ 
sented  in  this  figure,  the  temperatures  predicted  by  the  CMA  one-dimensional 
analysis  were  corrected  to  include  two-dimensional  heat  conduction  effects. 

Figures  4-20  through  4-23  show  the  measured  temperatures,  the  temperatures 
predicted  by  the  one-dimensional  analysis,  and  the  estimated  two-dimensional 
analysis  temperatures  at  the  backside  of  the  AGSR  for  nozzles  no.  4  through  7. 

These  comparisons  show  that  two-dimensional  conduction  effects  are  fairly  well 
characterized  by  the  temperature  correction  correlation  presented  in  Figure  4-19. 
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FIGURE  4-19  ONE-OIMENSIQNAL  ANALYSIS  TEMPERATURE  CORRECTION  INCREMENT 


Furthermore,  from  the  comparisons  shown  in  Figures  4-17,  4-18,  and  4-20 
through  4-23,  it  is  concluded  that  all  of  the  firing  and  cooldown  combinations 
could  be  accurately  predicted  by  the  2-D  analysis  but  that  useful  response  infor¬ 
mation  was  obtained  from  the  1-D  analyses.  The  primary  disadvantage  of  the  one¬ 
dimensional  analysis  is  that  it  underpredicts  the  backside  temperature  and,  thus, 
it  is  non-conservative  in  inferring  bond  line  temperatures.  In  addition,  the 
amount  of  discrepancy  between  the  temperatures  predicted  by  the  two  methods 
dictate  that  a  two-dimensional  analysis  is  required  if  a  valid  structural 
analysis  is  to  be  performed. 

4. 3. 1.2  Evaluation  of  the  Pyrolytic  Graphite  Layer  Recession  Calculations 

As  described  in  Section  4.2,  the  pyrolytic  and  ATJ  graphite  ablation 
rates  during  the  multiple  pulse  firings  were  evaluated  using  the  kinetically 
controlled  graphite  ablation  model  developed  under  Contract  F04611-69-C-0081. 

A  summary  of  the  computed  surface  response  of  the  pyrolytic  graphite  coated 
throat  insert  is  given  in  Table  4-9  for  each  firing  pulse.  Also  shown  in 
this  table  are  the  average  nozzle  throat  plane  surface  recession  measurements, 
average  throat  recession  rate  values  inferred  from  pressure  data,  and  the  pred¬ 
icted  temperatures  at  the  nozzle  surface  and  at  the  front  and  rear  faces  of 
the  AG SR  substrate.  For  nozzle  no.  8  both  the  one-  and  two-dimensional  heat 
conduction  results  are  shown.  As  pointed  out  in  Section  4. 3. 1.1,  the  two-dim¬ 
ensional  solution  gave  higher  temperatures  for  the  backwall  of  the  AGSR  subst¬ 
rate.  However,  the  surface  temperature  and  recession  results  were  only  slightly 
affected  by  the  2-D  heat  flow. 

The  average  surface  recession  comparisons  given  in  Table  4-9  indicate  that 
for  all  five  multiple  pulse  nozzles  the  predicted  throat  recession  was  substan¬ 
tially  below  the  average  of  the  post-test  recession  measurements.  For  nozzles  4, 
6,  7,  and  8,  the  PG  throat  liner  was  completely  burned  through  in  several  places. 
Only  the  nozzle  no.  5  throat  liner  which  was  exposed  to  only  two  pulses  survived. 
The  possible  reasons  for  the  difference  between  the  predicted  and  measured  surface 
recessions  are 

e  The  effect  of  nonuniform  combustion  product  mixing  on  the  local 
surface  chemical  reactions 

•  The  possibility  of  pyrolytic  graphite  surface  reaction  kinetics  chang¬ 
ing  because  of  the  annealing  effect  of  the  first  pulse  heating  cycle 

•  The  kinetically  controlled  reaction  rates  and  model  may  be  inaccurate 

•  The  difficulty  in  evaluating  the  heat  transfer  coefficient  values  in 
the  boundary  layer  separation  and  reattachment  region 

Aspects  of  each  of  these  possible  explanations  are  described  below. 
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*MTHMA  2-0  ANALYSIS 

As  described  in  Section  3,  there  was  evidence  that  the  combustion  products 
were  not  totally  mixed  (Figure  3-2) .  If  the  fraction  of  oxygen  in  the  boundary 
layer  were  on  the  average  higher  than  that  in  the  total  flow,  measured  recessions 
would  be  greater  than  the  computed  recession.  The  possibility  that  the  kinetic 
response  model  developed  under  Contract  F04611-69-C-0081  did  not  adequately  char¬ 
acterize  the  ablation  phenomena  in  the  multiple  pulse  duty  cycles  also  exists. 

For  example,  the  effect  of  the  pyrolytic  graphite  annealing  during  the  first 
firing  pulse  and  cooldown  may  have  affected  the  kinetic  response  mechanisms  during 
later  pulses.  Arc  heater  tests  were  run  under  the  kinetic  model  development  pro¬ 
gram  (Contract  F04611-69-C-0081)  which  indicated  that  annealed  models  tended  to 
ablate  slightly  faster  than  Has  deposited”  models.  Such  conditioning  of  the  PG 
coating  would  result  in  higher  recession  rates  during  the  start  pulses.  The 
recession  rates  inferred  from  chamber  pressure  versus  time  data  indicate  that  an 
annealing  effect  was  probably  not  important.  Although  the  quality  of  these 
recession  rate  data  is  low  because  the  changes  in  pressure  were  small,  the  average 
throat  recession  rates  were  roughly  constant  until  the  probable  burn  through  time, 
independent  of  the  type  of  start  up  or  the  number  of  pulses.  Alternately,  it  may 
be  possible  that  the  kinetic  response  model  is  not  suited  to  the  ANB  3066  simulator 
propellant  environment.  This  is  not  likely,  since  the  model  has  been  shown  to 
adequately  characterize  the  response  of  pyrolytic  graphite  to  several  similar 
propellant  formulations. 

The  most  probable  explanation  for  the  large  measured  throat  recessions  is 
the  severe  heating  caused  by  the  disturbed  flow  field  in  the  throat  region.  Boundary 
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layer  reattacbment  in  the  vicinity  of  the  throat  plane  could  provide  sufficient 
additional  heat  and  mass  transfer  to  explain  the  larger  than  computed  recessions. 

In  the  multiple  pulse  nozzles  the  shape  of  the  ablative  part  at  the  nozzle  entrance 
most  certainly  caused  a  separated  turbulent  wake  flow  field  in  the  throat  region. 
Reattachment  of  such  a  flow  would  result  in  a  very  severe  local  heating  condition. 
The  analytical  evaluation  of  the  heat  and  mass  transfer  in  such  a  flow  situation 
was  not  attempted  in  this  study.  However,  the  simple  flow  model  which  was  used 
for  the  computations  described  here  gave  a  baseline,  lower  limit  boundary  layer 
solution.  The  throat  package  thermal  response  comparisons  described  in  Section 

4. 3. 1.1  demonstrated  that  the  overall  boundary  layer  solution  was  adequate. 

However,  local  perturbations  such  as  increased  heating  at  the  reattachment  point 
would  not  significantly  affect  the  in-depth  response  of  the  nozzle  throat  insert. 

4.3.2  Carbon  Cloth  Phenolic  Exit  Cone  Response  Analysis 

One-dimensional  transient  heating  and  material  ablation  solutions  were 
performed  using  the  CMA  code  for  the  upstream  (A/A*  ■  4.2)  thermocoupled  plug 
location  for  each  of  the  five  multiple  pulse  nozzles.  The  1-D  analysis  procedure 
accurately  modeled  the  ablative  response.  Two  CMA  program  options  were  alter¬ 
nately  utilized  to  specify  the  heated  surface  boundary  condition.  For  each 
simulator  firing  pulse,  the  radiative  and  convective  surface  heat  fluxes  were 
evaluated  by  demanding  energy  conservation  at  the  ablating  surface.  These  solu¬ 
tions  accounted  for  the  energy  fluxes  associated  with  each  of  the  ablation 
phenomena  depicted  in  Figure  4-1.  For  the  cooldown  periods  between  or  after 
firing  pulses,  the  convective  heating  option  was  not  allowed  and  the  surface 
radiative  boundary  condition  was  evaluated  by  a  simplified  surface  energy  balance. 

The  nodal  system  for  the  CMA  analysis  was  the  same  for  each  of  the  five 
nozzles.  This  system  is  shown  schematically  in  Figure  4-24.  Also  shown  are  the 
locations  of  the  three  in-depth  thermocouples.  The  computed  temperatures  at 
these  locations  are  compared  to  the  response  measurements  in  Section  4. 3. 2.1. 
Recession  and  char  thickness  calculations  are  compared  to  post  test  measurements 
in  Section  4. 3. 2. 2.  Computed  density  profiles  are  related  to  the  density  pro¬ 
file  measurements  (Section  3)  in  Section  4. 3. 2. 3. 

4. 3. 2.1  In-Depth  Temperature  Histories 

The  computed  surface  and  in-depth  temperature  histories  for  the  five 
multiple  pulse  nozzles  are  compared  to  thermocoupled  plug  data  in  Figures  4-25 
through  4-29  The  general  conclusion  reached  from  these  comparisons  is  that  the 
analysis  technique  does  a  reasonable  job  computing  the  in-deptn  thermal  response 
of  the  ablative  materials.  Differences  between  the  response  analysis  results 
and  the  data  are  primarily  attributable  to: 
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•  uncertainties  in  material  thermophysical  properties 

•  the  inability  to  model  two-dimensional  heat  conduction  or  surface 
heat  convection  during  cooldown 

Factors  which  are  of  secondary  importance  in  explaining  these  differences  are 

•  uncertainties  in  thermocouple  depths 

•  erroneous  data  due  to  thermocouple  lag  time 

•  effects  of  char  swell  on  material  response 

An  analysis  of  the  sensitivity  of  computed  in-depth  temperatures  to 
material  property  errors  was  described  Section  4.2  which  showed  that  the 
property  values  used  in  the  1-D  analyse  were  of  sufficient  accuracy  to  compute 
the  proper  response  trends.  However,  uncertainties  in  material  properties 
would  result  in  an  uncertainty  band  of  at  least  +  100*F  in  the  predicted  temperatures. 

Differences  between  the  predicted  and  measured  temperatures  indicate  that 
the  material  response  was  not  well  modeled  during  the  cooldown  phases  of  the 
duty  cycles.  This  was  apparent  since  calculated  temperatures  in  the  char  layer 


4-41 


from  TC  no.  1  were  well  below  the  corresponding  thermocouple  data.  There  are 
two  possible  reasons  for  the  computed  results  being  low.  One  is  that  the  surface 
boundary  condition  specification  during  the  cooldown  was  incorrect.  However , 
this  explanation  does  not  seem  appropriate  as  the  analysis  of  the  temperature 
response  sensitivity  to  the  radiation  interchange  model  (Section  4.2)  indicated 
that  the  underprediction  could  not  be  explained  by  an  error  in  the  radiation 
model.  This  analysis  showed  that  even  with  no  net  radiation  energy  loss  from 
the  surface,  the  measured  temperatures  at  the  TC  no.  1  location  were  above  the 
computed  temperatures.  Furthermore,  it  is  possible  but  unlikely  that  during 
the  nozzle  cooldowns  the  simulator  purge  flow  provided  convective  heating  to  the 
nozzle  surface.  The  second  and  most  reasonable  explanation  for  the  measured 
temperatures  being  higher  them  the  1-0  solution  is  that  heat  is  being  conducted 
into  the  exit  cone  from  the  graphite  throat  package.  This  speculation  is 
confirmed  by  the  response  comparisons  for  nozzles  no.  4  and  5  (Figure  4-25  and  4-26). 
Each  of  these  nozzles  had  cold  restarts  and  the  difference  between  the  predicted 
and  measured  cooldown  temperatures  at  the  first  thermocouple  location  (0.125  inches) 
was  much  smaller  in  the  first  pulse  than  the  second.  This  £3  consistent  with 
the  realization  that  after  the  second  pulse  the  char  layer  would  be  much  thicker 
which  would  result  in  the  net  heat  flux  from  the  throat  insert  being  much  larger. 

The  heat  conduction  in  the  exit  cone  char  layer  is  difficult  to  model 
analytically  with  the  1-D  analysis  procedure.  A  crude  approximation  could  be 
made  by  allowing  a  source  term  to  be  specified  for  each  node  of  the  analysis 
grid.  A  more  rigorous  technique  would  be  to  use  the  1-D  charring  solution  for 
the  firing  pulse  as  the  initial  condition  for  a  2-D  non-charring  analysis  of  the 
cooldown.  The  2-D  grid  would  include  the  entire  throat  package  and  exit  cone. 

The  correct  modeling  of  cooldown  events  is  important  only  because  of  the  hot 
restart  duty  cycle  requirement.  For  hot  restarts,  the  errors  in  the  temperature 
field  are  maintained  and  subsequent  errors  become  additive.  For  the  multiple 
pulse  duty  cycles  analyzed  in  this  program,  errors  in  computed  in-depth  temperature 
were  no  larger  than  400*F. 

Tho  secondary  factors  described  above  may  also  have  contributed  to  the 
differences  between  the  computed  in-depth  temperature  histories  and  the  data. 

These  sources  of  error  are  not  peculiar  to  the  multiple  pulse  duty  cycles,  and 
since  they  have  been  discussed  in  detail  in  Reference  8,  only  a  brief  summary 
of  them  is  included  here. 

The  uncertainty  in  thermocouple  depth  directly  results  in  an  uncertainty 
in  the  predicted  temperatures.  Since  the  thermocouple  depths  reported  by  the 
nozzle  fabricator  for  the  thermocoupled  plugs  were  not  obtained  from  post  fabrica¬ 
tion  location  measurements  (e.g.,  from  x-ray  photographs),  the  uncertainty  in 
thermocouple  location  was  estimated  to  be  t  0.02  inches. 
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This  uncertainty  results  in  temperature  differences  of  approximately  +  200*F. 
Another  source  of  uncertainty  between  the  measured  and  predicted  temperatures  is 
the  thermocouple  la?  time.  This  lag  often  occurs  during  times  of  rapid  tempera¬ 
ture  rise,  particularly  if  tt;.*  thermocouple  junction  is  not  in  good  contact 
with  the  relative  material.  The  ovher  source  of  uncertainty  is  char  swell 
which  distorts  the  thermocouple  depths  end  makes  direct  comparisons  difficult. 
More  ie  said  about  the  char  swell  phenomenon  in  Section  4. 3. 2. 3. 

Despite  the  minor  difficulties  associated  with  the  effects  of  material 
property  uncertainties,  2-D  heat  conduction  effects  during  cooldowns  and  inherent 
systematic  measurement  uncertainties,  the  predicted  thermal  responses  of  the 
thermocoupled  plug  at  A/A4  *4.2  agreed  favorably  with  the  tempera ture  data. 

The  principal  conclusion  to  be  drawn  from  the  comparisons  between  the  predicted 
and  measured  in-depth  temperature  histories  was  that  the  thermal  response  of 
an  ablative  material  subjected  to  any  multiple  pulse  duty  cycle  is  well  charac  * 
terised  by  current  analysis  techniques. 


4. 3. 2. 2  Recession  and  Char  Thickness  Comparisons 

In  Table  4-10  computed  surface  recession  end  char  thickness  results  are 
compared  to  average  measurements  from  the  upstream  (A/A*  ■  4.2)  thermocoupled  plug 
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X-  ition  of  the  five  multiple  pulse  nozzles.  Because  the  material  was  a  carbon 
cloth  phenolic  fabricated  at  a  0  degree  layup  angle  with  respect  to  the  nozzle 
centerline/  char  swell  phenomena  had  to  be  considered.  Char  swell  apparently 
occur?  because  the  pressure  which  is  created  during  the  resin  depolymerization 
forces  the  distance  between  the  cloth  tape  lamina  to  increase  in  the  char  layer. 
These  forces  are  the  most  probable  cause  of  the  local  delamination  which  occurred 
at  several  places  in  the  exit  cone.  The  appearances  and  typical  locations  of  these 
delaminatione  were  indicated  in  Figures  3-1  through  3-5.  The  following  sketch  shows 
how  the  char  layer  thickness  is  increased.  Since  this  phenomenon  is  not  modeled  by 


the  CMA  analysis  procedure,  computed  results  and  data  cannot  be  conveniently  re¬ 
lated  without  a  char  swell  correction.  In  Reference  2  an  empirical  relation  was 
defined  so  that  the  predicted  thermochemical  surface  recession  values  could  be 
corrected  for  char  swell.  This  char  swell  correction  relation  modeled  the  obser¬ 
vation  that  the  char  swell  was  directly  related  to  the  char  thickness.  The  char 
swell  was  given  by: 


6 swell  *  keTchar 

Tchar  “  the  predicted  char  thickness  in  the  absence  of  char  swell 

kB  -  is  a  proportionality  factor  which  is  a  function  of  the 
a  function  of  the  material  and  the  local  layup  angle 

In  Reference  2,  the  proportionality  factor  was  found  to  be  0.135  for  a  NX  4926 
carbon  phenolic  exit  cone  part  with  a  layup  angle  relative  to  the  exposed  surface 
of  45  degrees.  The  0.135  factor  has  been  used  to  successfully  describe  the  char 
swell  amounts  for  a  variety  of  ablative  materials  and  layup  angles  (References  2, 
4,  and  5).  For  the  FN  5055A  carbon  cloth  phenolic  material  in  the  exit  cone  of 
the  multiple  pulse  nozzles,  the  char  thicknesses  were  computed  using  the  0.135 
char  swell  factor,  and  they  were  on  the  average  of  14  percent  lower  than  the 
measured  thicknesses.  Predicted  char  thickness  values  with  and  without  char  swell 
are  compared  to  data  in  Table  4-10.  The  computed  surface  recession  with  and 
without  the  char  swell  correction  are  also  compared  to  measurements  in  this  table. 


The  possible  reasons  for  the  inaccuracy  of  the  char  swell  correction  are 

thati 

•  The  computed  char  thickness  was  low  because  the  2-D  conduction  along 
the  char  layer  during  cooldowns  was  not  accounted  for. 

•  The  data  used  to  define  the  0.135  carbon  swell  factor  were  based  on 
char  swell  thickness  for  layup  angles  (with  respect  to  the  nozzle 
surface)  of  45  degrees  whereas  the  layup  angle  for  the  nozzle  exit 
cone  being  considered  here  is  15  degrees. 

The  first  item  above  is  the  most  significant  because  the  char  swell  correlation 
factor  was  defined  in  Reference  2  based  upon  a  predicted  char  layer  thickness. 
Presumably ,  that  baseline  charring  material  ablation  analysis  was  not  complicated 
by  2-0  heat  conduction  effects.  The  correlation  is,  therefore,  dependent  upon  an 
accurate  (no  swell  included)  prediction  of  char  thickness.  If  the  2-D  conduction 
effect  which  was  hypothesized  in  Section  4. 3. 2.1  had  been  taken  into  account  in 
the  analysis,  the  predicted  char  thickness  and  hence  the  computed  char  swell  would 
have  been  larger. 

The  other  item  in  the  above  list  is  a  possible  reason  why  the  char  swell 
correlation  developed  in  Reference  2  may  not  be  applicable  to  these  multiple 
pulse  duty  cycle  analyses.  The  major  shortcoming  of  the  correlation  is  that  it 
does  not  include  a  layup  angle  dependence  even  though  swell  primarily  occurs  in 
the  across- laminate  direction.  If  such  a  dependence  were  incorporated  in  the 
correlation,  it  is  probable  that  the  proportionality  factor  would  be  larger  than 
0.135  for  the  FM  5055A,  15  degree  layup  material. 

Despite  the  difficulties  in  evaluating  the  amount  of  char  swell,  the 
following  conclusions  about  the  analysis  techniques  and  the  material  survivability 
in  multiple  pulse  duty  cycle  applications  can  be  made:  1)  Predicted  surface 
recessions  from  a  design  point  of  view  are  conservative,  yet  realistic;  2)  char 
penetration  calculations  give  results  which  are  within  the  measurement  uncertainties. 

4. 3. 2. 3  Density  Profile  Comparisons 

The  measured  density  profiles  obtained  from  the  shave  and  weigh  and  x-ray 
transmission  techniques  were  presented  in  Figures  3-7  and  3-11,  respectively. 

Also  shown  in  Figure  3-7  are  the  predicted  density  profiles  which  were  obtained 
from  the  one-dimensional  CMA  predictions  presented  in  Section  4. 3. 2.1.  In 
addition,  the  predicted  density  profile  for  MX  4926  thermal  conductivities  is 
shown  in  Figure  3-7i. 
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All  the  comparisons  presented  in  Figure  3-7  between  the  predicted  and 
measured  profiles  show  that  the  thicknesses  of  the  fully  charred  and  of  the 
pyrolysis  zones  are  underpredicted.  However,  the  predicted  density  profile  when 
the  MX  492$  thermal  conductivities  were  used  in  the  analysis  is  in  better 
agreement  with  the  measured  profile  than  when  the  modified  FM  5055A  thermal 
conductivities  were  used.  This  difference  in  the  predicted  profiles  was  expected 
as  the  MX  4926  thermal  conductivites  are  higher  in  the  temperature  range  where 
the  majority  of  the  in-depth  material  degradation  is  occurring.  Even  though  the 
uncertainties  in  the  material  thermal  conductivity  in  the  pyrolysis  zone  may 
be  the  reason  for  the  discrepancy  between  the  predicted  and  measured  profiles, 
the  primary  reasons  are  believed  to  be 

•  axial  heat  conduction  in  the  char  layer  from  the  nozzle  insert  to 
the  exit  cone  analysis  location  during  nozzle  cooldown 

•  the  char  swell  phenomenon  which  was  evident  in  the  surface  recession 
and  char  thickness  measurements  presented  in  the  previous  section 

The  effect  of  *-he  axial  heat  conduction  would  be  to  increase  the  thickness  of 
the  fully  charred  zone,  and  the  effect  of  the  char  swell  phenomenon  would  be  to 
increase  the  thickness  of  the  pyrolysis  zone. 

The  measured  densities  in  the  fully  charred  and  pyrolysis  zones  using  the 
x-ray  transmission  technique  are  in  poor  agreement  with  those  obtained  using  the 
shave  and  weigh  technique.  For  example,  the  fully  charred  densities  from 
Figure  3-11  is  approximately  32  lbs/ft1  (*  0.5  gr/cm 3 )  as  compared  to  approximately 
70  lbs/ft1  for  those  presented  in  Figure  3-7.  However,  the  measured  distance 
from  the  surface  to  the  virgin  material  is  the  same  using  both  techniques  (0.45 
inches  from  Figure  3-7j  compared  to  0.47  inches  (12mm)  from  Figure  3-11). 

As  explained  in  Section  3.2,  the  porosity  of  the  sample  in  the  fully  charred  and 
pyrolysis  zones  is  the  primary  reason  for  the  discrepancy  in  the  measured  densi¬ 
ties.  Future  testing  using  the  x-ray  transmission  technique  should  try  to  mini¬ 
mize  the  effect  of  the  material  porosity  by  using  thicker  samples. 
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SECTION  5 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  principal  objectives  of  the  study  reported  here  were  as  follows: 

e  Obtain  and  report  ablative  response  data  for  materials  used  in 
multiple  pulse  duty  cycle  rocket  nozzles 

e  Verify  the  ability  of  current  analysis  techniques  tc  predict  the 
response  of  ablative  materials  in  multiple  pulse  duty  cycle 
applications 

Both  of  these  objectives  were  satisfied.  Specific  conclusions  formed  in  the 
course  of  the  program  are  as  follows: 

e  The  Philco-Ford  TMSO-EM  test  facility  successfully  simulated  the 
ANB-3066  propellant  within  a  few  percent,  although  flow  field  non¬ 
uniformities  appeared  to  contribute  to  asymmetrical  nozzle  ablation. 

•  During  each  firing  pulse,  chamber  pressure  data,  throat  insert 
backwall  temperature  data,  and  thermocoupled  plug  temperature 
data  were  obtained. 

•  Net  surface  recessions,  char  thicknesses,  density  profiles,  and 
X-ray  transmission  data  were  obtained  f^om  post-fire  measurements 
of  each  of  the  five  nozzles. 

•  Throat  package  backwall  temperatures  were  accurately  computed  using 
a  two-dimensional  thermal  analysis  procedure  (ARGEIBL,  ACE,  and 
ASTHMA  computer  codes) . 

•  One-dimensional  thermal  analysis  procedure  (CMA  code)  was  shown  to 
be  inadequate  for  evaluating  in-depth  temperatures  of  the  PG  coated 
throat  insert,  although  predicted  surface  temperature  and  recession 
rate  agreed  favorably  with  two-dimensional  results. 

•  PG  coating  recession  rate  was  substantially  underprcdicted  by  the 
kinetically  controlled  surface  response  model  developed  under  Air 
Force  Contract  F04611-69-C-0031 ;  the  underprediction  is  attributed 
to  the  inability  of  the  boundary  layer  analysis  procedure  (ARGF.ILL) 

to  predict  the  severe  heating  at  the  point  of  boundary  layer  reattach¬ 
ment  on  the  throat  insert. 
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•  One-dimensional  (CMA  code)  in-depth  temperature  predictions  for  the 
FM  S0S5A  carbon  phenolic  material  at  an  exit  cone  area  ratio  of  4.2 
were  in  general  agreement  with  the  thermocoupled  plug  data. 

•  Significant  heat  conduction  from  the  graphite  throat  package  along 
the  carbon  phenolic  char  layer  occurred  between  firing  pulses. 

•  Prediction  of  carbon  phenolic  recession  at  the  exit  cone  area  ratio 
of  4.2  was  confounded  by  the  char  swell  phenomenon,  but  agreement 
with  data  was  obtained .using  a  semi-empirical  char  swell  correction 
relation . 

•  In  light  of  the  char  swell  phenomenon  and  uncertainties  in  measure¬ 
ment  techniques,  the  predicted  density  profiles  at  the  4.2  area  ratio 
agreed  favorably  with  both  the  shave  and  weigh  and  the  X-ray  trans¬ 
mission  measurements. 

•  X-ray  diffraction  measurements  indicated  the  anticipated  increase  in 
graphite  crystallization. 

The  general  conclusion  to  be  drawn  from  the  efforts  described  in  this 
report  is  that  current  analysis  procedures  are  suitable  for  the  evaluation  of 
material  response  during  multiple  pulse  duty  cycles.  No  significantly  differ¬ 
ent  thermal  events,  other  than  increased  significance  of  two-dimensional  heat 
conduction,  occur  during  multiple  pulse  firings.  It  can  be  concluded  that  ma¬ 
terial  response  phenomena  are  not  strongly  dependent  on  the  type  of  duty  cycle. 

Specific  recommendations  relative  to  these  multiple  pulse  duty  cycle 
analyses  are  as  follows: 

•  Post-fabrication  X-ray  photographs  of  thermocoupled  plugs  should  be 
taken  to  accurately  define  thermocouple  location  relative  to  the 
exposed  surface. 

•  Utilization  of  the  separation  producing  lip  in  the  throat  inlet 
region  should  be  evaluated;  cold  flow  tests  and/or  boundary  layer 
analyses  should  be  performed  to  estimate  the  surface  boundary  con¬ 
ditions  at  the  throat  insert  reattachment  point. 

•  The  evaluation  of  in-depth  temperature  distributions  within  graphite 
throat  inserts,  such  as  are  necessary  for  thermostructural  analyses, 
requires  the  consideration  of  two-dimensional  conduction  events. 

•  Consideration  should  be  given  to  axial  heat  conduction  at  all  loca¬ 
tions  in  nozzles  exposed  to  multiple  pulse  duty  cycles. 

•  The  purely  theoretical  evaluation  of  surface  recession  for  exit  cone  char¬ 
ring  ablator  materials  requires  the  development  of  a  comprehensive  char 
swell  analysis  capability. 
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APPENDIX  A 


MEASURED  CHAMBER  PRESSURE  HISTORIES 

The  five  instrumented  nozzles  were  tested  by  Aeronutronic  in  their 
solid  rocket  motor  simulator  under  Air  Force  Contract  F04611-69-C-0039.  At 
the  conclusion  of  this  testing,  the  nozzles  and  the  reduced  data  were  delivered 
to  Aerotherm  for  post-test  analysis.  The  reduced  data  for  the  chamber  pressure 
histories  which  were  measured  during  the  five  multiple  pulse  data  cycle  tests 
are  presented  in  this  appendix  in  Figures  A-l  through  A-5.  The  test  date  is  also 
shown  on  each  figure.  In  some  figures  {for  example,  Figure  A-lb) ,  two  curves 
are  shown  for  the  ignition  transient.  The  solid  curve  represents  the  chamber 
pressure  in  the  o2  -  H2  start  motor  and  was  reported  because  the  pressure 
port  in  the  chamber  experienced  plugging  as  a  result  of  ablation  of  the  silica 
phenolic  liner.  The  dashed  line  represents  the  estimated  chamber  pressure 
during  ignition.  As  shown  by  this  dashed  line  and  by  the  solid  line  in  Figure 
A- la  which  was  the  measured  chamber  pressure  (no  plugging  occurred) ,  the  chamber 
pressure  ignition  transient  is  a  smoothly  varying  function.  The  oscillations 
recorded  by  the  transducer  in  the  start  motor  are  a  result  of  the  transient  flow 
rates  which  are  caused  by  closing  the  values  to  the  start  motor  and  by  opening 
the  values  in  the  fuel,  oxygen,  and  nitrogen  lines  to  the  combustion  chamber. 
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APPENDIX  B 


MEASURED  TRANSIENT  IN-DEPTH  TEMPERATURE  HISTORIES 

As  described  in  Section  2.2  and  in  Appendix  A,  the  five  thermocoupled 
nozzles  were  tested  in  the  Aeronutronic  .•■olid  rocket  motor  simulator  and  the 
reduced  data  delivered  to  Aerotherm  at  the  conclusion  of  these  tests.  The 
measured  temperatures  for  these  tests  at  the  locations  shown  in  Figure  2-1  are 
presented  in  this  appendix  in  Figures  B-l  through  B-5.  In  these  figures,  the 
curves  are  identified  by  Tgybgcript*  T^e  subscript  "A"  refers  to  temperatures 
at  the  backside  of  the  AGSR  graphite  and  the  subscript  "E"  refers  to  temperatures 
at  the  backside  of  the  ATJ  graphite  (refer  Figure  2-1) .  The  number  subscripts 
denote  the  thermocoupled  plug  number  which  are  correlated  with  the  nozzle  number 
using  Table  B-l.  The  number  following  the  dash  (i.e.,  the  number  one  in  T5_^) 
refers  to  the  thermocouple  number  in  each  plug.  The  number  one  refers  to 
the  thermocouple  nearest  the  inside  surface  of  the  nozzle  while  the  number 
three  refers  to  the  thermocouple  farthest  from  the  surface. 

TABLE  B-l 

DEFINITION  OF  THERMOCOUPLED  PLUG 
NUMBERING  SYSTEM  AND  THEIR  INITIAL  LOCATION 


Nozzle  No. 

Thermocoupled 
Plug  No. 

Initial  Area 
Ratio  Location 

4 

5 

4.29 

1 

6.66 

5 

7 

4.24 

3 

6.63 

6 

709 

4.26 
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APPENDIX  C 


DESCRIPTION  OF  AEROTHERM  COMPUTER  PROGRAMS 

The  following  sections  summarize  the  capabilities  of  the  Aerotherm 
computer  programs  which  were  used  in  making  the  analytical  predictions  for 
this  contract. 

C.l  THERMOCHEMICAL  STATE  CODE  (ACE) 

C . 1 . 1  General  Capabilities 

The  ACE  program  solves  for  the  equilibrium  or  kinetically  controlled 
chemical  composition  for  a  variety  of  open  or  closed  systems  of  arbitrary  chem¬ 
ical  composition.  The  ACE  program  is  an  outgrowth  of  the  chemistry  subroutines 
of  the  BLIMP  program.  The  equations  which  are  solved,  the  solution  procedure 
which  is  utilized,  and  the  program  characteristics  are  described  in  Reference 
C-l  and  are  summarized  in  References  C-4  and  C-ll. 

The  ACE  program  has  the  following  major  options: 

A.  Closed  System  Chemistry 

•  Evaluation  of  equilibrium  (or  kinetically  controlled  time  dependent) 
chemical  state  for  assigned  pressure,  elemental  composition  of 

up  to  three  component  mixtures,  and  either  enthalpy,  entropy, 
or  temperature.  Provision  is  made  for  readily  carrying  out 
isentropic  expansion  and  compressions. 

•  Solution  of  normal  and  oblique  shock  relations  to  provide  the 
state  of  the  gases  downstream  of  the  shock  and  the  isentropic 
stagnation  state. 

B.  Open  System  Chemistry 

•  Calculation  of  surface  mass  balances  to  determine  a  relationship 
between  normalized  char  recession  rate,  normalized  pyrolysis  gas 
rate,  surface  temperature,  and  pressure  while  considering  either 
equilibrium  between  the  char  and  gases  adjacent  to  it  or  while 
considering  selected  rate-controlled  surface  reactions. 

All  of  these  options  are  formulated  for  completely  general  chemical 
systems.  Consideration  of  any  molecular,  atomic,  ionic,  or  condensed  species 
requires  only  the  inclusion  of  the  basic  thermodynamic  data  appropriate  for 
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that  species.  These  data  are  obtained,  for  example,  from  the  JANAF  Thermochemical 
Data  Tables  and  include  curve  fit  constants  for  entropy,  specific  heat,  and  heat 
of  formation. 

All  options  generate  a  complete  state  of  the  system,  including  molecular 
composition  and  thermodynamic  and  transport  properties.  The  surface  state  op¬ 
tion  provides  additional  information  as  discussed  below. 

The  surface  state  option  of  the  ACE  program  contains  a  number  of  fea¬ 
tures  which  make  it  very  powerful  in  the  analysis  of  ablation  data  and  thus 
enables  the  determination  of  the  governing  surface  physicochemical  phenomena. 

In  the  first  place,  one  does  not  have  to  choose  a  priori  the  surface  species. 

For  example,  in  the  case  of  the  Apollo  material,  even  if  minor  constituents 
are  neglected,  the  char  surface  could  be  Si02*,  C*,  CiC*,  Si*  or  SijN^*. 

The  thermochemical ly  controlling  species  will  be  discovered  by  ACE  as  part  of 
the  solution  process.  Secondly,  it  is  possible  to  isolate  species  or  component 
gas  mixtures  from  the  system  or  to  consider  rate- control led  surface  reactions 
or  surface-catalyzed  homogeneous  reactions.  Thirdly,  each  condensed  species 
can  be  assigned  a  fail  temperature  above  which  it  cannot  serve  as  the  surface. 

This  latter  capability  can  be  used  to  represent  mechanical  removal  of  a  species 
or  removal  of  a  species  above  its  melt  temperature.  A  fourth  major  capability 
of  the  ACE  program  is  that  it  permits  consideration  of  unequal  diffusion  coef¬ 
ficients  as  well  r.s  unequal  heat  and  mass-transfer  coefficients  through  an 
application  of  the  same  approximation  for  binary  diffusion  coefficients  utilized 
in  the  BLIMP  program.  A  principal  limitation  of  the  ACE  program  is  that  the 
surface  is  considered  to  be  homogeneous.  To  illustrate,  the  surface  of  the 
Apollo  material  would  be  predicted  to  be  one  of  the  species  listed  above,  but 
not  a  mixture  of  two  or  more  of  these  species. 

The  surface  state  options  of  the  ACE  program  provide  a  char  recession 
rate  normalized  by  a  mass-transfer  coefficient,  and  other  information  needed 
to  perform  an  energy  balance  on  the  surface  of  a  charring  ablating  material, 
as  a  function  of  pyrolysis  gas  rate  normalized  by  the  same  mass-transfer  coef¬ 
ficient,  surface  temperature,  and  pressure.  It  thus  does  not  by  itself  con¬ 
stitute  an  ablation  prediction  tool.  In  the  first  place,  it  is  necessary  to 
specify  the  mass-transfer  coefficient  a. d  this  cannot  be  done  precisely  with¬ 
out  solving  the  boundary- layer  equations.  Secondly,  the  determination  of 
surface  temperature  requires  the  soLution  of  a  surface  energy  balance.  One 
procedure  which  is  used  for  ablation  predictions  is  to  generate  surface  state 
solutions  with  the  ACE  program  in  the  form  of  punched  card  output.  This  is 
then  used  as  input  to  the  CMA  program  (Option  1'  which  performs  the  surface 
enerqy  balance;  the  operation  of  this  coupled  abtation  prediction  is  discussed 
in  Section  C. 3.1. 
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C.1.2  Examples  of  ACE  Program  Use 


The  ACE  program  receives  daily  use  at  Aerotherm,  particularly  for  ablation 
computations.  Almost  all  of  the  references  to  be  cited  during  the  discussions 
of  the  ablation  codes  also  contain  descriptions  of  uses  of  ACE.  References  C-3, 
C-4,  C-6,  C-8,  C-9 ,  C-10,  C-ll,  C-12,  and  C-13-16  are  of  particular  interest  in 
this  regard.  Reference  C-17  describes  a  materials  screening  task  carried  out 
entirely  by  ACE  and  covering  a  very  wide  range  of  materials  and  environments. 

C.2  AEROTHERM  REAL  GAS  ENERGY  INTEGRAL  BOUNDARY  LAYER 

(ARGEIBL)  COMPUTER  PROGRAM 

The  ARGEIBL  program  calculates  bulk  transfer  parameters  for  non-ablating 
laminar  and/or  turbulent  boundary  layers  for  axisymmetric  or  flat  plate  geomet¬ 
ries  using  an  energy  integral  method.  The  essence  of  this  method  is  the 
assumption  that  the  flat  plate  relationship  between  the  local  boundary  layer 
energy  thickness,  0,  and  the  Stanton  number,  CH,  is  valid  also  for  flows  with 
a  streamwise  pressure  gradient.  This  assumption  permits  solution  of  the 
pressure  gradient  boundary  layer  energy  integral  equation  for  0,  and  then 
employed  again,  permits  calculation  of  the  pressure  gradient  Stanton  number,  CH, 
from  0.  Appendix  A  of  Ref.  c-2  describes  the  theoretical  background  of 
ARGEIBL  in  some  detail. 

The  ARGEIBL  program  can  treat  any  real  gas  chemical  system,  provides  for 
transition  from  laminar  to  turbulent  flow,  and  permits  axial  variations  in  wall 
temperature  and  free  stream  properties.  The  thermodynamic  and  transport  properties 
required  for  input  to  the  program  are  generally  obtained  using  the  ACE  program 
and  the  axial  pressure  gradient  from  a  flow  field  analysis.  Other  required 
inputs  include  (1)  boundary  layer  edge  thermodynamic  state  over  the  surface, 
and  (2)  geometrical  specification  of  the  vehicle  or  body  shape. 

Output  from  the  code  includes  boundary  layer  heat-transfer  coefficient, 
and  boundary  layer  thickness  parameters  (displacement  and  momentum)  as  a 
function  of  body  location.  The  heat  transfer  coefficient  distribution  is  a 
required  input  parameter  to  the  heat  conduction  codes  discussed  subsequently. 
Boundary  layer  displacement  thicxness  enables  evaluation,  in  an  iterative  manner, 
of  boundary  layer  edge  entropy  by  considering  the  shock  angle  through  which  the 
"edge  gas"  has  passed.* 

Reference  C-3  has  some  descriptions  of  typical  uses  of  the  ARGEIBL 
code  in  ablation  predictions. 


Experience  with  this  technique  at  Aerotherm  has  shown  near  absolute  convergence 
with  only  1  iteration,  i.e.,  a  first  guess  and  1  more  solution. 
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C.3  CHARRING  MATERIAL  ABLATION  (CMA)  PROGRAM 

C . 3 . 1  General  Capabilities  and  Solution  Procedure 

The  CMA  program  is  a  coded  procedure  for  calculating  the  in-depth 
thermal  response  of  a  charring,  ablating  material.  The  basic  physics 
included  correspond  to  simple  charring  of  the  form 

plastic  *  char  +  gas 

according  to  a  three  component  Arrhenius  rate  law  where 

p  -  r<pA  ♦  pb>  +  u-n  Pc  (c-i) 

and 


Usually  components  A  and  B  represent  constituents  of  the  resin  or  plastic 
binder  and  F  is  the  volume  fraction  of  resir.  in  the  composite  material,  while 
component  C  represents  a  more  refractory  reinforcement.  The  program  is  often 
used  for  graphite  and  other  non-charring  materials;  for  such  problems  the  de¬ 
composition  reactions  (C-2)  are  simply  by-passed  and  the  material  treated  as  a 
single  component  (r  =  0) . 

The  program  is  an  implicit,  finite-difference  computational  procedure 
for  computing  the  one-dimensional  transient  transport  of  thermal  energy  in  a 
three-dimensional  isotropic  material  which  can  ablate  from  a  front  surface  and 
which  can  decompose  in  depth.  Decomposition  (pyrolysis)  reactions  are  based 
on  a  three-component  model.  The  program  permits  up  to  eight  different  backup 
materials  of  arbitrary  thickness.  The  back  wall  of  the  composite  material 
may  transfer  energy  by  convection  and  radiation. 

The  ablating  surface  boundary  condition  may  take  one  of  three  forms: 

OPTION  1  -  Film  coefficient  model  convection-radiation  heating 
with  coupled  mass  transfer,  including  the  effects  of 
unequal  heat  and  mass  transfer  coefficients  (non-unity 
Lewis  Number)  and  unequal  mass  diffusion  coefficients. 
Surface  thermochemistry  computations  need  not  presume 
chemical  equilibrium  at  the  surface,  as  described  in 
detail  in  Section  c.1.1. 
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OPTION  2 


-  Specified  surface  temperature  and  surface  recession  rate. 

OPTION  3  -  Specified  radiation  view  factor  and  incident  radiation 
flux,  as  a  function  of  time,  for  a  stationary  surface. 

Any  combination  of  the  first  three  options  may  be  used  for  a  single 
computation.  Option  3  is  appropriate  to  cooldown  after  termination  of  con¬ 
vective  heat  input  and  is  often  useful  in  conjunction  with  Options  1  and  2. 

The  program  permits  the  specification  of  a  number  of  geometries,  in¬ 
cluding  plane,  cylindrical  or  annular,  and  spherical.  In  the  most  general 
case,  area  may  vary  arbitrarilv  with  depth. 

The  rear  surface  of  the  last  node  may  be  specified  as  insulated,  or 
may  experience  convective  and  radiative  heat  transfer  to  a  "reservoir"  at  a 
specified  reservoir  temperature  if  a  rear  surface  convective  coefficient  and 
an  emissivity  are  input. 

Material  properties  such  as  thermal  conductivity,  specific  heat,  and 
emissivity  are  input  as  functions  of  temperature  for  virgin  plastic  and  char. 

For  partially  decomposed  material,  the  program  performs  an  appropriate  averag¬ 
ing  on  density  to  determine  effective  material  properties. 

The  basic  solution  procedure  is  by  a  finite  difference  approach.  For 
each  time  step,  the  decomposition  relations  are  solved  and  then  the  in-depth 
energy  fluxes  constructed  in  general  terms.  These  are  then  harmonized  with  a 
surface  energy  balance  (  if  a  surface  energy  balance  option  is  being  used)  and 
the  in-depth  temperatures  determined.  New  material  property  values  are  set  up 
and  the  solution  is  ready  for  the  next  time  increment. 

The  CMA  program  outputs  instantaneous  mass  ablation  rates  and  blowing 
parameters  for  char  and  pyrolysis  gas,  total  integrated  mass  ablation  of  char 
and  pyrolysis  '■'is,  total  recession  and  recession  rates  of  surface,  of  the  char 
line,  and  of  the  pyrolysis  line.  It  also  outputs  the  surface  energy  flux  terms, 
namely,  the  energy  convected  in,  energy  radiated  in,  energy  reradiated  out, 
chemical  generation,  and  conduction  away  (<3con{j)  •  Further,  it  describes  how 
the  input  energy  of  <3cond  is  "accommodated"  or  "partitioned"  in  the  solid  mat¬ 
erial.  Part  of  the  energy  is  consumed  in  decomposing  the  plastic,  part  is 
consumed  in  sensible  enthalpy  changes  of  the  solid,  and  part  is  "picked  up” 
by  the  pyrolysis  gases  as  they  pass  through  the  char.  Thermocouple  and  iso¬ 
therm  output  can  also  be  called  for. 
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C. 3.2  Some  Surface  Energy  Balance  Details 
C. 3.2.1  General  Remarks 

The  most  powerful  and  unusual  feature  of  the  CMA  program  is  its  general 
coupling  to  a  thermochemical  erosion  energy  balance  at  the  exposed/  heated 
surface.  This  coupling  involves  two  aspects:  the  surface  energy  balance 
formulation  itself,  described  in  this  section,  and  the  chemical  state  pro-ram 
used  to  provide  input  information  for  the  energy  balance  process,  described  in 
Section  C.l. 


C. 3.2. 1.1  Surface  Energy  Balance  Procedure 

In  calculations  under  Option  1,  the  in-depth  solution  is  coupled  to  a 
general  film-coefficient  boundary  condition. 

The  following  sketch  shows  the  various  energy  fluxes  of  interest 
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Since  the  surface  control  volume  is  attached  to  the  receding  surface,  there  is 
an  apparent  flow  of  char  material  into  the  control  volume  and  an  associated 

In  the  case  of  a  charring  material,  the  pyrolysis  gas  flow 

There  is  a 


energy  flux  mchc 


into  the  control  volume  from  below  bringing  an  energy  flux  m  h  . 
diffusional  energy  flux  in  from  the  boundary  layer;  this  is  denoted  qdiff 
There  is  a  convective  energy  flux  away  from  the  surface  in  the  gross  motion  due 
to  mass  injection:  (pv)whw.  An  additional  important  gross  motion  energy  flux 
is  denoted  £m  *h,  in  the  sketch;  this  term  represents  energy  carried  away  in 
the  removal,  through  mechanical  action,  of  condensed  phase  materials. 

Other  terms  in  the  energy  balance  are  the  obvious  radiation  terms  and  conduction 
into  the  material. 


In  the  course  of  the  transient  in-depth  solution,  which  essentially 
provides  the  <lconcj»  it  is  necessary  to  evaluate  all  of  the  other  surface 
energy  balance  flux  terms.  This  in  turn  naturally  requires  some  sort  of 
surface  thermochemical  state  solution.  For  this  purpose  it  has  proved  ex¬ 
pedient  to  prepare  in  advance,  using  suitable  ablation  thermochemistry 
codes,  (such  as  the  one  described  in  Section  C.l)  a  scries  of  tables  which 
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include  all  the  surface  mass  transfer  and  chemical  relations.  The  in-depth 
solution  may  then  be  coupled  to  the  surface  events  through  the  surface  energy 
balance.  For  example,  when  chemical  equilibrium  is  achieved  at  the  ablating 
surface,  the  development  presented  in  Reference  C-4  describes  the  means  for 
obtaining  the  thermodynamic  state  of  the  gas  at  the  ablating  surface  in  terms 
of  the  pressure,  and  char  and  pyrolysis  off-gas  rates.  In  equation  form. 


Thermodynamic  state  =  f(B',B‘,P) 

9  c 


(C-3) 


where 

m 

B'  =  o~u*^C"  {normalized  pyrolysis  off-gas  rate) 

9  pe  e'M 

B'  =  n~"T ;V  (normalized  char  recession  rate) 

c  °e  e  M 

P  =  boundary  layer  edge  pressure 

The  thermodynamic  state  includes  definition  of  surface  temperature  and  gas 
molecular  composition.  This,  in  turn,  enables  evaluation  of  the  various  quan¬ 
tities  appearing  in  the  boundary  layer  driving  potential  for  heat  and  mass  trans 
fer  (Reference  C-4) .  Tables  representing  solutions  to  the  functional  relation¬ 
ship  represented  by  Equation  (A-3)  are  generated  for  a  complete  map  covering 
the  range  of  B^,  B^,  and  P,  of  interest.  These  tables  are  generated  by  the  ACE 
chemical  state  program  described  in  Section  C.l. 

As  an  example  of  this  procedure,  suppose  a  table  is  prepared  which,  for 
a  parametric  array  of  dimensionless  char  erosion  rates  (B^) ,  dimensionless  pyrol 
ysis  off-gas  rates  (B^) ,  ana  pressure,  presents  the  relevant  ablating  surface 
temperature  and  requisite  boundary  layer  composition  and  enthalpy  quantities. 
During  each  time  step  in  the  course  of  the  in-depth  solution  the  program  gener¬ 
ates  a  pyrolysis  gas  rate  B^  and  computes  the  rate  at  which  energy  is  conducted 
into  the  material  from  the  surface.  Kith  and  the  pressure  known,  the  input 
parametric  tables  then  serve  to  define  that  which  yields  temperature  and  en¬ 
thalpy  quantities  which  provide  a  balanced,  harmonized  set  of  energy  fluxes  at 
the  surface. 

It  is  worth  noting  that  the  coupling  between  one  in-depth  solution 
(performed  by  the  CMA  program)  and  the  surface  thermochemistry  solution 
(performed  by  the  ACE  program)  yields  a  very  general  and  powerful  ablation 
thermochemical  data  to  bo  extended  to  many  geometries  and  flow  fields  at 
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no  additional  cost,  since  the  surface  thermochemical  solution  quantities  are 
normalized  on  the  transfer  coefficient.  Secondly  and  more  importantly,  the 
ACE  program  is  general  enough  in  its  thermochemical  solution  procedures  to 
handle  any  combination  of  environment  and  ablator.  Thus,  the  CMA  plus  ACE 
scheme  can  treat  complex  chemical  systems  as  easily  (from  our  standpoint)  as 
-the  simple  carbon-air  system,  in  practice,  these  pregrams  have  been  used  with 
great  success  as  heatshield  and  rocket  nozzle  design  and  analysis  tools  for  a 
wide  variety  of  materials  (including  graphite,  carbon  phenolic,  graphite  phen¬ 
olic,  silica  phenolic,  tungsten,  tantalum,  ice,  and  teflon)  exposed  to  varied 
environments  (including  aluminized  solid  propellants;  beryllium  and  solid 
propellants;  liquid  propellants  such  as  H2~C>2 ,N204-N2H4/UDMH,  OF2-B2Hg,  and 
f lox-methane ;  various  hybrid  propellants;  as  well  as  air  and  numerous  simple 
test  gases  such  as  Ar,  He  and  N2-  Furthermore,  ACE  will  automatically  discover 
the  chemical  nature  of  the  exposed  surface  (  it  may  be  a  non-obvious  species 
containing  elements  from  the  environment  as  well  as  from  the  ablating  materia1), 
and  can  account  for  any  number  of  kinetically  controlled  heterogeneous  reactions. 
These  two  features  are  very  important  ones  for  graphite  ablation  analysis. 
Graphites  are  sometimes  "loaded"  with  "additives"  to  form  erosion  resistant 
coatings  on  the  heated  surface  (as  described,  for  example,  in  Reference  A-5) . 

ACE  will  predict  the  existence  of  these  coatings  and  the  resultant  changes  in 
erosion  resistance.  Kinetic  control  is  occasionally  of  interest  for  graphite 
also.  ACE  will  also  include  the  effects  of  unequal  mass  diffusion  coefficients 
for  the  various  boundary  layer  species. 

ACE  can  alto  automatically  include  simple  melting  of  surface  species, 
detecting  surface  species  and  allowing  them  to  run  off  at  their  melting  points, 
accounting  for  the  energy  events  correctly  in  the  process  (Emr*h^  in  the  sketch 
preceding) .  Thus  one  form  cf  thermomechanical  ablation  (that  occurring 
as  a  total  loss  of  strength  at  a  given  temperature  for  a  given  species)  is 
automatically  accounted  for  in  CMA  through  the  ACE  solution.  Other  forms  of 
mechanical  removal  can  be  directly  accounted  for  in  CMA.  This  task  is  particu¬ 
larly  simple  for  a  non-charring  material  such  as  graphite  since  the  mechanical 
failure  rate  does  not  in  fact  change  the  surface  energy  balance;  the  failing 
mass  rate  both  enters  and  leaves  the  surface  energy  balance  control  volume  at 
the  same  temperature.  Thus  the  thermochemical  corrosion  rate  can  be  determined 
in  the  usual  manner  and  the  surface  temperature  calculated  in  the  process. 
Following  this,  the  thermomechanical  erosion  rate  can  be  determined  using  what¬ 
ever  procedure  seems  best  (such  as  with  empirical  correlations)  and  the  total 
recession  rate  determined  as  the  sum  of  thermochemical  and  thermomechanical 
events . 
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C.3.3  Uses  of  the  CMA  Program 

The  CMA  program  has  received  extensive  use  in  the  Option  2  (assigned 
surface  temperature  and  recession  rate)  mode  to  match  in-depth  thermocouple 
data  and,  in  doing  so,  to  discover  in-depth  properties  such  as  thermal  conduc¬ 
tivity.  References  C-6  and  C-7  provide  examples  of  such  use.  The  CMA  program 
has  also  been  very  extensively  used  in  its  Option  1  mode  as  an  ablation  pre¬ 
dictor;  References  C-3,  C-6,  C-8,  C-9,  C-10  provide  a  few  descriptions  of  the 
many  uses  of  the  program. 

C . 4  THE  AXI-SYMMETRIX  TRANSIENT  HEATING  AND  MATERIAL  ABLATION 

ABLATION  PROGRAM  (ASTHMA) 

C.4.1  Introduction 

For  strongly  two  dimensional  ablation  problems  such  as  nose  tips  and 
certain  rocket  nozzles,  1-D  calculations  of  the  kind  performed  by  the  CMA  pro¬ 
gram  described  above  must  be  supplemented  by  two  space  dimension  ablation  calcu¬ 
lations.  The  ASTHMA  program  is  a  suitable  2-D  code  for  this  purpose.  It 
accounts  only  for  non-charring  materials  such  as  graphite;  otherwise  it  is  as 
general  in  all  respects  as  the  CMA  program  described  above.  In  particular,  the 
-ode  has  the  same  very  general  thermochemical  erosion  boundary  condition  with 
appropriately  general  space  dependent  input  of  transfer  coefficient,  pressure, 
and  boundary  layer  edge  or  recovery  enthalpy  so  that  2-D  thermochemical  erosion 
may  be  realistically  predicted.  The  following  sections  describe  the  ASTHMA  code. 

<_' .  4 . 2  General  Capabilities 

The  ASTHMA  program  computes  the  transient  temperature  response  within  a 
two-dimensional  axisymmetric  non-charring  material.  It  accounts  for  two-dimen¬ 
sional  surface  recession,  and  allows  for  anisotropic  in-depth  thermal  conduction. 
Input  consists  of  geometry  specification  (in  terms  of  finite  difference  grid), 
material  properties  tables,  and  heated  surface  boundary  condition  information 
as  a  function  of  time  and  surface  location.  The  ASTHMA  program  has  the  same 
three  boundary  condition  options  as  the  CMA  code  described  above: 

OPTION  1  -  Film  coefficient  model  convection-radiation  heating  with 

coupled  mass  transfer,  including  the  effects  of  unequal  heat 
and  mass  transfer  coefficients  (non-unity  Lewis  number)  and 
unequal  mass  diffusion  coefficients.  Surface  thermochemistry 
computations  need  not  presume  chemical  equilibrium  at  the 
surface . 

OPTION  2  -  Specified  surface  temperature  and  surface  recession  rate. 
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OPTION  3  -  Specified  radiation  view  factor  and  incident  radiation  flux, 
as  functions  of  time,  for  a  stationary  surface. 

Any  combination  of  the  first  three  options  may  be  used  during  the  problem  history 
of  any  surface  point. 

Program  output  includes  the  history  of  the  in-depth  temperature  field, 
the  history  of  the  surface  point  locations  (recession  and  shape  history)  and 
surface  temperatures,  and  the  history  of  the  local  surface  energy  and  mass 
fluxes . 

C.4.3  The  Surface  Energy  Balance 

As  in  the  case  of  the  CMA  program,  the  most  powerful  and  unusual  feature 
of  ASTHMA  is  its  general  coupling  to  a  thermochemical  erosion  energy  balance 
at  the  heated  surface.  This  energy  balance  is  essentially  the  same  as  used  in 
the  CMA  code;  it  is  described  in  detail  in  Sections  C.3.2  above.  Here  it  will 
sufficient  toemphasize  that: 

1.  This  procedure  will  treat  any  material  exposed  to  any  environment. 

2.  It  accounts  for  unequal  diffusion  coefficients,  non-unity  Lewis 
number,  departures  from  chemical  equilibrium,  and  mechanical  ablation. 

3.  The  procedure  will  discover  any  chemical  surface  coatings  naturally 
formed;  this  is  important,  for  example,  for  loaded  graphites. 

4.  Boundary  condition  input  (recovery  enthalpy,  transfer  coefficient, 
pressure,  and  radiation  flux)  may  be  made  a  function  of  surface 
location. 

C. 4 . 4  U3es  of  the  ASTHMA  Code 

The  ASTHMA  code  is  described  in  References  C-13  and  C-I4,  where  graphite 
rocket  nozzle  throat  insert  ablation  studies  are  briefly  described.  Reference 
C-15  gives  more  detail  on  these  studies.  Reference  016  describes  numerous 
applications  of  ASTHMA  to  graphite  nosetip  ablation  predictions. 
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Measurements  and  analysis  of  the  thermal  and  ablative  response  of  five 
multiple  pulse  duty  cycle  rocket  nozzles  were  performed.  Nozzles  were 
fired  on  a  solid  propellant  simulator  with  a  test  stream  composition, 
designated  ANB-3066,  which  contained  16  percent  aluminum.  Nominal  cham¬ 
ber  pressure  for  each  of  the  various  firing  pulses  was  750  psia.  Materic 
temperature  histories  were  recorded  at  several  locations  in  each  nozzle. 
Post- test  measurements  of  surface  recession,  char  penetration,  and  den¬ 
sity  versus  depth,  as  well  as  thermocouple  data  are  used  for  comparison 
with  one  and  two-dimensiona]  heat  conduction  and  material  ablation  calcu¬ 
lations.  Agreement  between  prediction  and  data  is  favorable,  although 
the  analyses  showed  that  two-dimensional  heat  conduction  effects  need  be 
accounted  for  to  obtain  an  accurate  evaluation  of  the  thermal  response 
of  multiple  pulse  duty  cycle  rocket  nozzles.  It  is  concluded  that  mate¬ 
rial  decomposition  and  ablation  phenomena  are  not  strongly  dependent  on 
the  type  of  duty  cycle. 
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